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1*0 INTRODUCTION 


This interim report covers the activities of this Grant NSG-4027 for the 
third quarter of the grant year which ends 31 May 1983. This third quarter 
covers the period from 1 December 1982 through 28 February 1983. The 
following personnel are assigned to the grant research working 1 /4-time under 

l 

the grant funding: 


Dr. Richard K. Smyth 
Mr. Phillip Chan 
Mr. Fadi J. Kurdahi 
Mr. David Ho 
Mr . Carposf oro Sosa 
Mr. Jean-Francois Soi 


Principal Investigator 
Research Assistant 
Research Assistant 
Research Assistant 
Research Assistant 
Research Assistant 


Effective 1 June 82 
Effective 1 Jan 83 
Effective 1 June 82 
Effective 1 Jan 83 
Effective 1 June 82 
Effective 1 Jan 83 


The assignments of the research assistants are as follows: 


(1) Design, Code, & Test Air Data System Software and Investigate 
Hardware; Carposf oro Sosa and Fadi Kurdahi 

(2) Complete Hardware & Software for Spin Warning System Designed 
by Team #1 Utilizing SC-01 for Voice Generation; Philip Chan 

(3) Complete Software, and Test the On-Board Simulation; Jean- 
Francios Soulard 


(4) Provide Support in Software Coding, Compilation, and Object 
Code Loading into 68000 Microcomputer for all the Grant Research 
Teams; David Ho 


In addition to the personnel assigned to the grant, other graduate 
students taking the EE560L microcomputer research course, and students taking 
directed research EE590L under Dr. R.K. Smyth, have performed research which 
have contributed to the grant's technical objectives. The final reports of 
these student team's research reports are attached as appendices to this 
Quarterly Interim Status Report. The Graduate Students and their contribu- 
tions to the grant objectives follow: 


Team #1 (Report Appendix A) 


1. On-Board Simulator 

Mr. Jeffrey Bluen 
Mr. Jean-Francios Soulard 
Mr. Mehdi Namakian 
Mr. Charles Saleh 

2. Spin Warning System Team #1* (Report Appendix B) 

Mr. David Barry *used SC-011 voice generation chip 

Mr . David Ho 

Mr . Renshan Tang 

Mr . Mohammed Movahed-Ezazi 

3. Spin Warning System 

Mr. Steve Meier 
Mr. Tieh Ku 
Mr. Tom Wilkenson 
Mr . Dave Adachi 
Mr. David Chen 


Team #2** (Report due May 83) 

**used TI voice generation chip 



4. On-Board Simulation 

Mr. Jeffrey Bluen (EE590L) 
Mr. Horng-Ru Hwang (EE590L) 


Team #2 (Report due May 83) 
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2.0 AIR DATA SYSTEM (ADS) 


The ADS equations defined in section 2.4.1 of the Semi-Annual Interim 
Status Report (30 November 1982) have now been coded, and the cooding for the 
equations are being tested. 

The work on ADS covers two sub-areas: installation of software support 
tools and implementation of the air data system software. With respect to the 
first sub-area, during this period we have finished the installation and 
testing of the software tools resident on the IBM 370/4341 (ECL-VIRGIL) . This 
means that written software for the three areas of research (ADS, SWS, and 
OBS) , can now be tested in the VERSAM0DULE-01 system. 

Upon completing testing of the M68000 cross-software (Jan 25) we proceeded 
to continue our work in the air data system. Currently we are in the area of 
developing software to test each of the 13 implemented ADS equations in order 
to obtain an estimate of their execution time. These estimates are required 
to design an appropriate time scheduling. 

The progress in the two sub-areas is explained in detail in the following 
sections . 

2.1 Software Support Tools 

The system software support tools are applicable to all of the tasks on 
the grant, although the ADS team is checking out the tools. The following 
tasks have been completed. 

TASK 1 : The debugging of the communication software (CMSCPM) that handles file 
transfers between IBM (VIRGIL) and Coraputerm (EE560L) has been completed. A 
short user's guide, describing the use of CMSCPM was written. 

TASK 2 : Modification of the linker so that it produces a load map of the 
execution module being generated. A version of this map, displaying the 
absolute addresses of all the sections comprising the module, is required to 
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identify the entry point, i.e. the address at which execution starts. This 
task was achieved by modifying the EXEC (CMS command file) that invokes the 
linker, inserting new options and the corresponding file definitions. 

TASK 3 : Test the execution of PASCAL programs in the VERSAMODULE-Ol . In this 
respect we encountered two problems: 

(a.) Incorrect loading of the execution module from Computerm into 
VERSAMODULE-Ol. The second line of the S-format file, containing 
the execution module, did not get loaded into the monoboard 
memory. This file, resident in Computerm, is transferred from the 
IBM using the utility CMSCPM, already mentioned in Task 1. We 
compared the files in the two systems and we did not detect any 
modification caused by transmission errors. We also checked the 
header of the S-format file and replaced it with others, used for 
S-format files that are known to load normally. The fault 
persisted. We have avoided this problem by duplicating the second 
line using a text editor. 

(b.) Identification of the section in the execution module that 
contains the entry point. Motorola supplied us with the correct 
entry point for a module running under EXORMACS but did not give 
any information for running a program under VERSAMODULE-Ol. With 
the aid of the road map, we surveyed each section until we found 
the correct entry point. 

We have now checked the test file supplied by Motorola as well as some 
preliminary programs used in ADS. So far, we have not encountered any 
problems regarding their execution in the VERSAMODULE-Ol. 

2.2 Air Data System Software 

The 13 ADS equations (table 2-1) have been coded. The ADS equation timing 
tests (table 2-2) have been completed and the timing procedure software has 
been written. The data acquisition routines have been completed and the 


testing of this software is underway. A potentiometer test fixture is being 
designed to permit testing of the various ADS algorithms which will emulate 
the sensor voltage levels with the pots and test for correct computation of 
the various air data parameters such as Mach number, altitude, etc. 


I ’ 
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TABLE 2-1 

AIR DATA SYSTEM PARAMETER EQUATIOHS 

(Page 1 of 3) 


PTI (TOTAL PRESSURE) : 
PTI = QCI + PSI 


MI (INDICATED MACH NUMBER): 
for PTI/PSI <= 1.893 

MI = SQRT (5.0) * SQRT( ( (PTI/PSI) ** (2/7)) -1) 
for PTI/PSI > 1.893 

X = 1.839371 * (PSI/PTI) 

MI = SQRT ( (Ax - Bx - Cx* *2 - Dx**3 - Ex**4 - Fx**5 - Gx**6 - 
Gx**9 ) /x) 


MINF (FREESTREAM MACH NUMBER) 

MINF = MI + DM 

DM is an error correction obtained by interpolating 
look-up table 


in a 



PSINF (FREESTREAM STATIC PRESSURE) 
for MINF <=1 

PSINF = PTI / ( (1 + 0.2 * MINF* *2 ) ** (7/2)) 
for MINF >1 

PSINF = (PTI * A * (l-A) ** (5/2) / 0.1839371 
where A = 1 / (7 * MINF**2) 


QBAR (DYNAMIC PRESSURE) 

QBAR = 0.7 * MINF**2 * PSINF) 

QCC (CORRECTED AIRSPEED PRESSURE) 
QCC = PTI - PSINF 


KEAS (KNOTS EQUALIVANT AIRSPEED) 

KEAS = MINF * 661.48 * SQRT(PSINF / 2116.22) 


KCAS (KNOTS CALIBRATED AIRSPEED) 

KCAS = 1479.1 * SQRT ((1 + (QCC / 2116.22) ** 


(2/7)) -1) 


(- 


* . origWal pase is 
r: of ipoor quality 

TABLE 2-1 

AIR DATA SYSTEM PARAMETER EQUATIONS 

(Page 2 of 3) 

HP (GEOPOTENTIAL OR PRESSURE ALTITUDE) 
let R = PSINF / 2116.22 
for R > .223361 

HP = 145442 * (R ** .1092632365 - 1) 
for .223361 => R > .0540328 

HP = 164219.39 - 20805.7 * Ln(PSI) 

for .0540328 => R > .00856663 

HP = 710793.96 * A**2 - 645177.17 
where A = (.0540328 / R) ** .01463563358 
for R <= .00856663 

HP = 81660.714 * A* *2 - 162928.85 

where A = (.00856663 / R) ** .04097977402 


AINF (ANGLE OF ATTACK) 


f (MINF) 
EA 

MINF 

.0055 

V 

o 

. 

o 

.0053 

0.2 

.0051 

0.4 

.0044 

0.6 

.0033 

o 

« 

CO 

.0023 

0.9 

.0 

V 

j— ■ 
. 

o 


AINFF = (1 + EA) * (ALPHA I) 

BINF (ANGLE OF SIDESLIP) 

BINF = (1 + EB) * (BETA I) 
where EP = 0.0 


, original page is 
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TABLE 2-1 

AIR DATA SYSTEM PARAMETER EQUATIOHS 

(Page 3 of 3) 


GAMMA (FLIGHT PATH ANGLE) 

GAMMA = THETA - A INF 

HDGAMMA (ALTITUDE RATE) 

HDGAMMA = 60 * MINF * CS * Sin (GAMMA) 

where CS = Speed of sound per 1962 std. atmoshpere 

HDOT (ALTITUDE RATE / TIME DERIVATIVE 
let HAV(t) - 1/5 HP(t-i); 

HDOT(t) = HAV(t) - HAV(t-l) [time interval = 1] 

FQTY (FUEL QUANTITY) 

FQTY (0) = 304.85 AND FUSED (0) = 0.0 

if FFR(t) => 3.25 and FFR(t) <= 52.0 then FUSED ( t ) = FFR ( ti 
if FFR(t) < 3.75 or FFR ( t)' > 52.0 then FUSED (t) = FUSED (!• 

and 

FQTY (t) = FQTY (t-1) - FUSED (t) 

when TOPOFF = 1 FQTY(t) is reset to 304.85 

once LAUNCH = 1 never reset FQTY(t) 
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TABLE 2-2 

TIMING TESTS FOR ADS COMPUTATIONS 
(See Table 2-1 for ADS Parameter Equations) 


Execution time in milliseconds 


PARAMETER 

PARAMETER NAME 

TEST 1 

TEST 2 

TEST 3 

AVERAGE 

PSINF * 

Free Stream Static Pressure 

85.70 

81.33 

87.48 

84.84 

QBAR 

Dynamic Pressure 

5.115 

5.115 

5.115 

5.115 

QCC 

Corrected Airspeed Pressure 

2.13 

2.10 

2.07 

2.10 

KEAS 

Knots Equivalent Airspeed 






( true ) 

17.58 

17.23 

17.17 

17.33 

KCAS 

Knots Calibrated Airspeed 






(indicated) 

88.03 

90.33 

86.71 

88.36 

HP ** 

Geopotential or Pressure 






Altitude 

55.95 

170.38 

87.12 

104.48 

AINF 

Angle of Attack 

14.145 

14.145 

14.145 

14.145 

GAMMA 

Flight Path Angle 

2.130 

2.10 

2.055 

2.10 

HDGAMMA 

Altitude Rate (Computed for 






GAMMA) 

39.43 

40.48 

40.48 

40.13 

FQTY 

Fuel Quantity 

3.165 

3.150 

3.150 

3.16 


TOTAL 

313.38 

426.36 

345.50 

361.75 


MAX ITERATION RATE 3.19/sec 2.35/sec 2.89/sec 2.76/sec 


SAFE INTERATION RATE 2/ sec 


* Two equations, selected by freestream Mach No. value 
** Four equations, selected by freestream Mach No. value 
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The student team cited in the introduction and lead by Bluer & Soulard 
produced a final report on the OBS which is included as Appendix A of this 
progress report. Messrs. Bluen and Soulard are continuing work on the OBS. 
Mr. Hwang is providing coding for the aircraft lateral equations of motion (3 
degree-of-f reedom) using the T-38 aircraft parameters for the equations. He 
is using the data from a NASA research report written by Mr. Teper of Systems 
Technology, Inc. (STI) for the equations and parameters. His program will be 
menu-driven and will permit the selection of various flight conditions to be 
used. Messrs. Bluen & Soulard will integrate Mr. Hwang's equations and 
software modules into the overall OBS software. 

The modification and improvements to OBS team # l's project are described 
below. The system concepts are covered in figures 3-1 through 3-5. 

3.1 Iapleaentatioa of Separate Motions 

In the earlier form of this project the line-of-sight and the range 
between the two planes were generated by programmed functions. This choice 
implied that the host and the target were not actually moving independently. 
The new implementation generates separate geometry parameters for both air- 
planes. These parameters are updated during every time step of the simulation 
loop. The program calculates the Cartesian coordinates of the two planes (XH 
and YH for the host and XT and YT for the target) using the velocities 
(VELOCITY and TVELOCITY respectively) and the turning angles (PSI and PSIT 
respectively) given by the simulation loop. The variations are first 
calculated (DELTA-XH, DELTA-YH, DELTA-XT, DELTA-YT) and then added to the old 
values of the coordinates (XOH, YOH, XOT, YOT). At the beginning of the loop 
the newly computed values are assigned to the old value variables and the 
process is started again for the new time step. 

r 

In addition, both aircraft use separate simulation loops to generate the 
turning angle commands. 
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3.2 Hew Tracking Procedure 

The "prediction" algorithm presented in the former report proved to be 
inappropriate for the use of the OBS. It has been replaced by a conventional 
guidance law, called the proportional guidance system, in which the parameters 
are calculated using actual geometry of the scene (as opposed to "anticipated" 
position of the target as used in the "prediction" algorithm). The parameters 
used by the proportional guidance are: 

■ The line of sight rate, SIGMA 

■ The guidance factor, LAMBDA 

■ The guidance gain, GUIDGN 

These parameters are calculated with the updated geometry given by the 
program, that is, passed to the guidance law that generates an optimal command 
for the host. The pursuit loop can be closed automatically (as is done now) 
or by an actual pilot using a display screen. 

3.3 Providing Target Maneuvers 

At the beginning of the simulation, the pilot or the operator is asked the 
initial geometry of the scene. He must give the coordinates of the target 
relative to his starting position. The host starting position is taken as the 
origin of the grid. The target maneuver is deterained by its initial turning 
angle command, or alternately, can be programmed as a sequence of commands. 

3.4 Providing Different Types of Aircraft 

Two possibilities are considered: 

■ Providing fixed pre-defined types, chosen in a menu, or 

■ Providing ad libitum types under reasonable limits 

Choosing pre-defined types could ease the initialization procedure for the 
user but, on the other hand, ad libitum types allow an infinite range of 
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FIGURE 3-4 ANALYTICAL BLOCK DIAGRAM FOR OBS 


H = HOST AIRCRAFT 
T = TARGET AIRCRAFT 
K t = ANGLE TRACKER GAIN 

n~ T = ANGLE TRACKER TIME CONSTANT 


OF POOR QUALITY 
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aircraft but require more work from the user. An operation manual may be 
suitable but a decision has not yet been made. 


3.5 Proportional Navigtion 


Proportional navigation tries to keep the proper lead angle from host to 
target aircraft always forcing the host line-of-sight rate, CT , to zero. It 
determines the host aircraft velocity vector turning rate which is propor- 
tional to the line-of-sight rate. The host aircraft velocity heading, , 
changes until R <T becomes zero. The most important component of this 
computation is the accurate measurement of the spatial line-of-sight rate. 


The aircraft is assumed to have a large accurate antenna used to give 
line-of-sight information and is assumed to be space stabilized using free 
gyros. We will further assume that the antenna is slaved to keep its axis 
perfectly aligned with the gyro axis, so we will not have to model gimbal and 
torquer dynamics. Initially, no noise on the steering signal is assumed. 


The proportional relationship between O'h and is constantly readjusted by 
the closing rate as determined from the changing geometry. The relative 
velocity is divided by the host aircraft velocity and multiplied by the 
assumed_A. of 4.0. 


X 

Vr 


*1 

a 



steering gain 
closing rate along R 
host aircraft closeing rate 
along R 
line-of-sight 
look angle 

velocity vector direction 
turning rate 


7 aos^-t 
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4.0 SPIN WARNING SYSTEM 


Mr, Philip Chan fixed Team #l's breadboard voice generator using the SC-01 
chip. During early February he had the breadboard speaking words and phrases 
using phonemes programmed into the system using a CRT keyboard. He is 
programming the works and phrases contained in the SWS matrix shown in table 
4-1. The phoneme method of generating words used by the SC-01 chip appear to 
be a flexible and feasible method. 


TABLE 4-1 


STALL WARNING & COMMAND TO RECOVERY 


Absence of 

Rate DC ±10 Discrete 1 Discrete 2 Discrete 3 Discrete 1/2,3 


YAW RATE 

ASYMMETRIC THRUST l 


AB BLOWOUT 

AB STALL 

MIL STALL 

SYMMETRIC THRUST 
(NO BLOWOUT OR STALL) 

Less than 
40°/sec 

"Unload" 

If IAS < 150 KTS 
"Left or Right 
(Lighted AB) 
Engine Mil" 

"Unload" 

"Left or Right 
(Stalled) Eng Idle" 
If IAS 4175 KTS 
"Left or Right 
(Unstalled) Enq Mil" 

"Unload" 

’ If IAS 4 150 & 
Altitude 4 15K 
"Left or Right Eng 
‘ Idle" (Unstalled Eng) 

No Warning 

Less than 
50°/sec but 
greater than 
40°/sec 

"Unload" 

"Left or Right 
Engine Idle" 

"Unload" 

"Both Engines Idle" 

"Unload" 

"Both Engines Idle" 

"Unload" 

Greater than 
50°/sec 

"Left or Right 
(Lighted AB) 
Engine Mil" 
"Stick Full Left 
or Right" 

"Stick Full Fwd" 

"Both Engines Idle" 
"Stick Full Left or 
Right" 

"Stick Full Fwd" 

"Both Engines Idle" 
"Stick Full Left or 
Right" 

"Stick Full Fwd" 

"Stick Full Left 
or Right" 

Push Stick 
into turn 


IAS 


After any warning from table above: 

If Yaw Rate<40°/sec for two second "Center Controls". 

If Yaw Rate <30°/sec for two seconds & Airspeed >120 Kts 
A0A<-10° and Ay> .1g "Stick Half-Aft and Hold": then; 

If A0A> 0 for 5 sec "Recovery Complete". 


DC 

±10 


0-10 


'Recovery Complete" 


Stick Pos Pitch ±1011 
Roll ±I0U 
Angle of Attack 
Lateral Acceleration 

±10V to 50K 
= TN g's 1g 
Indicated Airspeed 
(in knots) 


A0A = 

\ = 
ALT = 
UNLOAD 
IAS = 


o o 
■n 2 
*0 P 
O 2 

8P 
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Either engine above 1215° (EGT turbine limit temp) for 3 sec "Left" or "Right Engine Off" 
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5.0 ACTIVITIES PLANNED FOR NEXT (LAST) QUARTER 

The Air Data System software should be completed and tested using simul- 
ated sensor signals through the analog-to-digital converter card. The ADS 
software should be available for NASA Dryden to use on their flight simulator 
by the end of May 1983. 

The on-board simulator will integrate the three degree-of-freedom (3 DOF) 
equations with the intercept algorithms. Emphasis will be placed on 
generalizing the OBS to tasks other that the air-to-air intercept mission 
currently programmed. 

The Spin Warning System effort being done by the team #2 will provide a 
final research report to be incorporated in the Grant Final Report produced in 
May 1983. Mr. Chan will continue working on the Team #1 and will implement 
the PASCAL-based software on the 68000 microcomputer and integrate the 
software with the SC-01 voice generator. 


APPENDIX A 


ON-BOARD SIMULATOR 
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INTRODUCTION 


> t 


K'u:,n 


The project presented here is a synthesis of two 
subjects suitable for EE 560L, the graduate research course 
for Advanced Microcomputer Applications at the University of 
Southern California. The two original parts from which this 
work derives are: 

On Board Flight Simulator (Project #13.0) 

Air-to-Air Intercept Mode (Project #18.0) 

The On Board Simulator part contains a three degree- 
of-freedom Aircraft Behavior Simulation, providing parameters 
used by the Interception procedure. These parameters could 
also be used for verifying closed loop performance before 
flight . 

The Air-to-Air Intercept Mode is a software package 
integrated in the simulation process that generates a Target 
motion and performs a tracking procedure that predicts the 
most likely next target position, for a defined time step. 

This procedure also updates relative position parameters and 
gives adequate fire commands. The simulation of the input 
data, provided by an angle tracker, is done by pre-chosen 
programmed functions. This allows a wide range of target 
behavior as well as the full control of it when testing the 
procedures. 
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ON BOARD SIMULATOR 


Preventing mishaps and saving tremendous amounts 
of money and time are the motivators for the construction of 
this 16-bit microcomputer simulation. The immediate benefits 
are derived from the fact that the On Board Simulator (OBS) 
is intended to use the existing hardware - aircraft flight 
computer. It can be easily reprogrammed to simulate different 
aircraft, and it also uses an inexpensive processor. It can 
provide parameters to check out the aircraft flight computer 
as an added side benfit. 

The OBS in effect can use part of the host Aircraft's 
electronics to close the loop for the simulator thus spreading 
the computing load. A simulator of this type previously 
required the computing ability of much larger units and so 
large simulations were built justifying the greater computer 
expense. With the advent of the small, efficient, and computa- 
tionally powerful 16-bit microprocessor, a new approach to 
simulation is possible. It can now execute the requisite 
equations in a real time sense enabling flight simulation in 
this hybrid configuration. 

The On Board Simulator will be connected to the 
actual aircraft via the onboard computer and pilot commands. 

It will receive aircraft flight commands from the pilot, will 
simulate the aircraft dynamics and return to the aircraft 
flight computer updated parameters of the vehicle geometry. 
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Various aerodynamic and geometry equations can be programmed 
depending upon the type and model of aircraft being modeled. 

The proposal to EE 560L included a 3 DOF of 
aircraft dynamics which would receive information from the 
pilot joystick as he appraises the engagement geometry seen 
on his video display. His input was a command to the host 
aircraft roll rate. 

The approach to building the. On,, Board Simulator 
began with the high level language, PASCAL - , ; t?h'ree degree-of- 
freedom simulation of the dynamic parameters. The working 
model was designed and tested and results are included. The 
remaining step is to download to object code in the Motorola 
68000 microcomputer. This has not yet successfully been done 
because of unavailability of the required software and hard- 
ware tools. 

A preliminary program designed to test the CRT 
driver and A/D conversion is also included. It was not fully 
implemented. (See section on preliminary program.) 


AIR-TO-AIR INTERCEPTION 

The air-to-air interception procedure was implemented 

using different modules: a tracking process using both the 

parameters passed by the simulation and a target motion 

generator, that could be reprogrammed at one's will. According 

» 

to the data provided by both the generator and the simulation. 


the next target position of the target is estimated and the 
optimal next position of the host aircraft is calculated. 
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AIRCRAFT 


DYNAMICS 


7 , 
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The focus of this project is the simulation of both 
the Host (H) and the Target (T) planes motions. The Host 
plane's pilot gives orders to his machine and we had to 
recreate the response of the instruments and the structure 
to these commands. That is why we need some fundamental 
Aircraft Dynamics results. 

In this introducing short study, we only try to 
show very basic data about Aircraft maneuvers and the para- 
meters on which the pilot can act. The equations and their 
computing are treated in the project. 

1 • THE PILOT'S COMMANDS : 

To control the plane's motion, the pilot has three main 
mechanisms : 

The THROTTLE, commanding the thrust, that is 

the power given by the engine. For a definite 
aerodynamic configuration of the plane, at a 
constant altitude Z=Zq, the throttle commands the 
speed of the plane. 

The stick, commanding the pitch, attack and 

roll angles by its action on the wing flaps. 

The rudder bar, commanding the rudder angle 

that directly acts on the yaw angle. 



These three mechanisms interact very closely, in such 
a way that it is hard to point out the effect of one of 
them without looking over the two others, (see Figure 1) 


• V' rw 01 ^figure I : 

CONTROL UNITS OF A : SIMPLE AIEPLANE . 
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BASIC EQUATIONS FOR AIRCRAFT DYNAMICS - CONSTANT ALTITUDE: 




Our method ofi/^l^baini'ng'.coprdination is based on 
the fact that for a certain bank angle and true air speed, 
there is only one value of yaw rate (\jj) for which coordina- 
tion can be achieved (refer to objective of the project. 
Aircraft Dynamics, 2.2). 



C (center of the Circle) 

llocfl K 


V T = tangential aircraft velocity 


V ? 
T i 2 

F sm<J) = m — — 

K 


=m ^ 2 R=mV T ^ 


F t coscj) = m g 


V T • 

=> tg <p = — ip => 


i 

= S~ tg<p 
T 

Approximation used: 
for <j> < 20° tg <p - (p -*■ 
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STUDY OF TURNING: 
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When a pilot wants his plane to turn, at a constant 
altitude, he must roll his aircraft. In order to cancel the 
effect of gravity the centrifugal force should be as shown below: 


»> 


FI ( Ur'T) 


Lift projection compensating 



w = rotation speed 

r = ray of the circle (or part 
of circle described by 
the plane) 


BASIC EQUATIONS: 


F = F => 
Ly c 


F l sintj) = m w‘ 


1 

2 


p S V 2 C , with: 
2 


F Lz = mg = > 


F l cos 


lj) = mg 


> m cos|) > n z cos 


f 


p = volumic mass of air at Z=Z o 

S = lifting area of the plane 

V = speed of the plane 

Cz= aerodynamic coefficient of 
lift 
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A PRIORI LIMITATIONS TO fTHE^'SIMULATlON MODEL; 

Within the aerodynamic equations we can see that for 

a given roll angle, the turning rate is entirely defined. 

That is why our model will only consider the roll rate as a 

variable and our pilot will only use the stick to control his 

airplane. Note that this is true only because we assume that 

the host aircraft and target aircraft are co-altitude. 

We also have to be aware of the priori limitations 

existing on a turning manuever. The human body cannot stand 

more than a limited acceleration, measured in number of G's 
2 

(1G = 9.81 m/s ). This number is represented in the turning 
equations by n sometimes called the charge factor. In a 
turning configuration we showed that n was equal to 1/cosi|j; 
the following table gives the values of turning angles for 
different charge number, measured in units of G. 


Value of n (G's) 
z 

Value of the turning angle * 

8 

84 degrees 

7 

83 degrees 

6 

80.5 degrees 

5 

79.5 degrees 

4 

77 degrees 

3 

75.5 degrees 


* Note that these angles are lowered by Aerodynamic considera- 


tions not taken into account here. 
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BUILDING THE CONTROL LOOP 


The pilot acts on the stick by giving it an angular 
displacement Sg T , that drives the wings ailerons to an angular 
displacement 6 that creates a banking moment. This moment is 
proportional to the rolling acceleration $, and the diagram 
leading to <p is the following: 


PILOT! 


LU 


When we have 4> , we can use the formula: \p = ~ tg cj> , 

V T 

• • 

to calculate ip. For this, we need first to integrate (p then 

we assume that tg4> ~ <J> , which is true for the small angles, and 

our formula becomes: ip = <p . The diagram between (f> and iji is: 

V T 


1 

. 

! 4> 

i 

_J g 

s 

| 



J 

1 T . 


These 

diagrams 


to feedback loop 

These diagrams do not take into account the delays 
that can occur in the commands. We tried to figure it by using 
the following control diagram: 


PILOT 


T 3 S + 1 


with the pilot being modelized as: 


K i (T i S + 1) 
T Z S + 1 


E 




kernel ■+/"'\ 


Ki (T i S+l ) 


Hi! I 


— C5- 


This is the final control loop description for the 
three degree-of-f reedora simulation in PASCAL. It has an added 
rate damping path with gain and is shown to be stable in the 
following root locus analysis. The first order lag acting on 
the stick command models the time delay in the dynamic reaction 
of the aircraft to the commanded rolling moment. 

The lag was previously placed with the aerodynamics 
but concern has arisen over the aircraft roll angle during its 
maneuvers. For accurate modeling and evaluation of possible 
limiting it was determined that the lag should delay the roll 
rate and not the turning rate. 

In evaluating the roll limiter it is noted that with- 
out any limiter on the roll angle it grows to 4.2 radians at 
time equals 1.5 seconds. This is clearly unrealizeable so a 
shunt limiter was placed on roll angle limit. Attached plots 
of heading, roll angle, and roll rate for both unlimited and 
limited cases clearly depict this situation and solution. 



L 
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Starting with the system: 





In order to determine a suitable gain (K^i^) a root locus analysis was 
performed on the open loop transfer function KGH: 

K (S + 2.5) (S + .5) 

KoH — — _ 

S 2 (S + 4) (S + 2) 


Where Ti = .4 



K^/K 5 = .5 

Kn = 1 

V = 200 m/s 
g =10 m/s 2 


Two methods were used to determine the optimal gain, the first is a hand 
calculation and the second corroborating analysis was performed using an 
interactive computer aided design tool called TOTAL which derives from the 
USAF E_glin base. 


KGH = 


K (S + 2. 5) (S + .5) 
S 2 <S + 4) (S + 2) 


1. Find the root locus asymptotes. 


OA = 


4 + 2 - (2.5 + .5) 6-3 


1 
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2. Find the root locus "real" value as it crosses the line + j . 

The phase angle criterion yields: for (-1+j) 

- [tan * 
18.4 


[90° + tan " 1 (-^2.) ] 


+ [ tan 


-l 


1.5 


] - 2 [135°] - [45°] 


+ 116.6 


33.7 


270 


- 45 


- 183. | 0 


3. Use the magnitude condition to find K at (l±j) 


iKGHl => K = 


GH 


1 S*a-l+j 


= (/TTT ) 2 /~3 2 ''~+~~P r /~2 
/ 1 + (1.5) ^ / 1 + . 5 * 

- 2/~2 /To 

/ 3.25 / 1.25 
= 4.44 


4. To find the roots that correspond to K = 4.44 

Choose a value on the real axis and compute its K value. 

Say S = -3.5 

= (-3,5) 2 (-3.5 + 4) (-3.5 + 2) = 12.25 (.5) (+1.5) = 

(-3.5 + 2.5) (-3.5 + .5) (1) (3) 

This is too small an S value 


3.06 


W|H> 
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Say S = -3.0 


(-3 . 0) 2 (-3.0 + 4.0) (-3.0 + 2.0) = 9 » 1 » 1 

(-3.0 + 2.5) (-3.0 + .5) ” (.5) (2.5) 

This is too large an S value 


Say S = -3.35 


(3. 35) 2 (+3.35 - 4.0) (3.35 - 2.0) 
(3.35 - 2.5) (3.35 - .5) 


11.22 (.65) (1.35) „ 

" (’.85) “27 85) = 4 ‘° 6 


S = -3.3 


K = 


(3.3) (-3.3 + 4.0) (3.3 - 2.0) 

(3.3 - 2.5) (3.3 - .5) 


10.89 (.7) (1.3) 

(. 8 ) ( 2 . 8 ) 


= 4.42 


Trying a similar point for the other pole zero pair which moves an equal 
distance: 


S = -1.26 


(1 . 26) 2 (-1.26 + 4.0) (-1.26 + 2.0) 

(-1.26 + 2.5) (1.26 - .5) 


1.59 (2.74) (.74) 

(1.24) (.76) 


3.42 


S = -1.0 


13 1 


(1.5) (.5) 


4.0 


S = -.9 


.9 2 (3.1) (1.1) 

(1.6) (.4) 


4.3 


S = -.78 


■78 2 (-.78 - 4.0) (-.78 + 2.0) _ .608 (3.22) (1.22) 

(-.78 + 2.5) (-.78 + .5) = (1.72) (.28) 


4.96 




Since K 
G 


V T 3 T 2 
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and K = K 5 

n 


If poles are at -3.3, -.9, (lHHj) K = 4.44 

What is the gain for design? 

K 4 44 m — Li — * 5*^*2 g T 1 - - 4 
T2T3 V 

.4 2 K 1 K 2 

.25 .5 ZQ 
10 

= .04 • 8 KiK 2 = .32 K 1 K 2 


4.44 

.32 


KjK 2 


13.875 = 


KiK 2 
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Checking by plugging in values: yji >' »•- 


G = 


K 

A 

KiK 2 g 

V S 2 (T 2 S + 1) (T 3 S + 1) 


(T l S + 1 ) 


K Ti (S + Tl) 

TaT?(S + 1 ) (S + 1 ) 

t 2 t 3 


H = K S (S + 

*5 


4(S) _ G 
4 c (s) 1 + GH 


KT ( (S + — ) 

3j 

T 3 T 2 S 2 (S + ^-) (S + ^L-) + K Tl (S + Ks (S + £j) 


K T| (S + Ti) 

T3T2 S 2 (S 2 + 7- + ~) S + — ) + K K S Ti (S 2 + (r 1 - + |^) S + -— 3-7 
Tl T 2 T1T2 — __ - Ti K5 T iK 


T.T, S' + * -1-) K_K i _ Ll S 2 + X + X ♦ 


-> K - 13 • 9 X 10 - figs 

K 200 " - 695 


X = K • K 5 Ti = . 695 « (2)x . 4 = .556 
■ K -— *— (S + — ) 

. T 3— L2 1j 


S“ +(7- + 7-) S 3 + &- + K K S Tj) S 2 + K K S T, (“f~ + ~) S + K K 5 Tj 

TlT 2 Ti T1T5 T1K5 


(.695) (.4) _l 

(.5) (.25) .4 


J3 

s'* + (-*r + ^r-) s 3 + .4 

• J * « J 


+ .556 s2 + .556(-^ + .5) 


.556 (- ^ ) 
S + .4 


(.5) (.25) 


(.5) (.25) 


(.5) (.25) 


2.22 (S + tW) 

s“ + 6 S 3 + 14.4 S z + 13.3 S + 5.56 


This is similar to the computer's TOTAL answers 
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The user package called TOTAL is an easy way to analyze control loops. One 
enters in , then/f(s) and from there performs requested options interactively. 
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TOTAL • • VERSION 3.0 
STANLEY J. LARIMER 


******** COLIN USERS ******** 
QUESTIONS OR PROBLEMS » CONTACTS 
CART DENNIS DIDALEUSKY 
AD/SDES-A (2-5678/2-5669) 


The general options categories available are: 



S AN INTERACT COE COMPUTER-AIDED DESIGN PRO ORAM 
ITAL S CONTINUOUS CONTROL SYSTEM ANALYSIS. 
TAINS 160 OPTIONS DIVIDED INlO GROUPS OF 10 
DING TO GENERAL APPLICATION. 


*J.0N3 ENDING IN O LIST 
ft EXAMPLE f OPTION SO LI 


THE NEXT 10 OPTIONS* 

STS OPTIONS 30 I HRU 39. 


THE FOLLOWING ARE TUI MAIN OPTION GROUPS: 

0-91 TRANSFER FUNCTION INPUT OPTIONS 
10-19,* MATRIX INPUT OPTIONS 
20-29J BLOCK DIAGRAM MAN! PULAT l ON OPTIONS 
30-391 TIME RESPONSE OPTIONS 
40-49: ROOT LOCUS OPTIONS 

50-59: FREQUENCY RESPONSE OP I IONS 

60-69: POLYNOMIAL OPERATIONS 

70-79: MATRIX UPS PAT IONS 

30-09: DIGITIZATION OPTIONS 

90-99: OPTIONS OF PARTICULAR 

100-109: MURE TRANSFER FUNCTION 
110-119: MORE MATRIX OPTIONS 
.120-129: MORE BLOCK DIAGRAM MAN IPUI. A TIONS 
130-139: STATE TRANSITION' SIMUI ATI ON OPTIONS 
140-149J DOUBLE-PRECISION DISCRETE IRANSFOPM OSPTJf'NS 
150-159J MULT I -RATE FREQUENCY RESPONSE OPTIONS 


INTER*' 3 ( 
.INPUT OF 


T IONS 


Y 
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The input open loop transfer function is as follows: 


Oil-- (if i'K/OUiK ' • 

GAIN- 1.00 0 


P . 0 yO 


OLIF(S) NUMERAIuK 

!;l .NFTJLY < I > OLZERO(I) 

2 000 )Stl 2 < - . t;ooo ) + J< y. 

6.000 ' )S.N 1 ( -2.500 > + i.' 

2*500 ; 01 NK- . . vO 

OLTF(S) DENOMTNA fuR 

HI.DPOLY < I > OLPOLE ( I ) 

l .000 )S** 4 <0. > »• ,J< 0. 

6.000 >S** 3 (0. ) + J( 0. 

0.000 >S** 2 < -2.000 ) + J< 0. 

>, >S ** 1 ( -4.000 ) + J( 0. 

) ol m<« i.ooo 


and the root locus options are: 


OPT ION 
* stO 
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l 

l • 

4t 

* 

42 


43 

* 

44 

»/ 

'I* 

45 

1* 

1' 

46 

* 

47 

Is 

T- 

48 

*i* 

•A 

49 


i . ■ • ' r locus up r rows 

r ISAM 1 ' v POP S-! i nfltr TSAMP*SAMPLE TIME FOR Z-PLANE ) 
1.1 > i OPTION'S 
ui.Ni-.krtl- «<uuT LOCUS 

KUO f LOCUS WITH A CAIN OF INTEREST 

ROOT LOCUS WITH ZEIA (DAMPING R A f I 0 ’ OF INTEREST 

LIST N PUINIS ON A BRANCH OF INTEREST 

LIST ALL POINTS ON A BRANCH OF INTEREST 

LIST LOCUS KUOTS AT A OAIN OF INTEREST 

LIST LOCUS ROOTS AT A ZETA OF INTEREST 

PLOT ROUI LOCUS AT USER'S TERMINAL 

LIST CURRENT VALUE'S OF ALL ROD I LuCUS VARIABLES 


m 


TYPE': HELP * 49 FuR DEFINITION OF ROOT 1.0CU3 VARIABLES 

A CALCUHP PI OT FOR OPTIONS 41 » 42 f 43 >4 40 MAY PL 
OBTAINED BY TYPING: PI. 0T,41 OR PI. 07 » 42 ETC. 


OPTION 
? 41 


30 


! H 'EN-T.f)t »F-' ( DLTF ) ROOT LOCUS USING OPTION 


40..; n". 


W 


4 I'Gl.rG A"! 


X 

X 

X 


0 . 

0 . 

o.oooo 

-4.0000 


Y = 

Y =•■ 

Y - 

Y = 


■ < 5 <' 
' 
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0 . 

0 . 

0 . 

0 . 


2 ZL'fc os 
x =■• 
x == - 


AT 

. 50000 


4- • sJ 


000 


Y = 

Y = 


0 . 

0 . 


CAIN CONSTANT < OLNN/OLDK > “ 


2.0000000 


REGION or CALCULATION-REAL: CC* -5.00 TO AA-* 1.00 
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The four branches plotted. root locus are tabulated enabling choosing 

a desired gain or damping -ratio. 

nz&tori op calcnlat nu.' real: cc- tc> i..-.o 

W /AVU1 DO- 3.1-u TO $8= 3 . 00 

A [{RANCH UlARIiNG AT (0.1 1 i(0.) 

• TYKE L i:i LIST , S TO Shir. OR T> TO ABORT :• 



i. 

original page fs 

of POOR QUALITY 

. :!\V NCI I NUflPLR L 

TLCULATION L ,‘fP Cl EE = ,1000 

- MINTING STEP SIZE = .1.500 



l. f=r*!'? REAL 

LOCUS JrtAG 

BIST TO ORIGIN 

t .AIN 

ZETA 

r.n 

( *. 4 


O . 

. tooooooor 02 

0 . 
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♦ 
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4 mI.TO 7 

X » •*. ■' U ( J 
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, 

t 
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. -.SJ 40? 2 
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nnUNDARY 

BRANCH STAR TING AT ( 0 . ) + J < 0 . ) 

TYRE J. TO LIST* S TO SKIP, OR 't TO ABORT > 


MRANCH STARTING AT ( 7 . > K J<0.) 

TYPE L TO L iSI» S TO SKIP* OR * TO ABORT > 


BRANCH NUMBER 3 


CALCULATION STEP SIZE * .3000 

PRINTING STEP SIZE = .1.500 



LOCUG REAL 

LOCUS IMAG 

BIST TO ORIGIN 

GAIN 

ZETA 

C B 

2.0000000 

0. 

7,0000000 

0. 

1.00000 

1 

1 . 3000000 

0. 

1 .9000000 

-* r> -r **> 0 *■? 

« t / / 

1,00000 

0 

1 .6C0GC00 

0. 

1 .6000000 

1 . 2-1121 

1.00000 

0 

1.4000000 

0. 

1 , 4000000 

1 .51 421 

1.00000 

O 

1.2000000 

0, 

1 . 2000000 

4 T'V'l'TI 

X ♦ / / «rUl 

3.00000 

0 

1 .0000000 

0. 

1 .0000000 

2.00000 

1,00000 

•0 

.30000900 

0. 

.90900000 

2.-10741 

1.00000 


.60000000 

0. 

.60000000 

4,50947 

3 ,00000 

Ti 

• SOO'oOOOO 

0. 

,50000000 

0. 

1.00000 



BRANCH STARTING AT (-4.) + J(0.) 

TYPE I. ID LIST TU SKIP* OR % W ABORT > 
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DRAMS! 1 

. j* in i.< KR 4 



CAL LU*. A T ION '•) i < r" ST. '».! - 

. 10 00 



PRINTING 

SILL SIZE 

. 1500 



l.OC'JS REAL 

LOCUS I MAG . 

HIST TO ORIGIN 

GAIN 

ZETA 

i .0000000 

0. 

. 0000000 

0. 

1.00000 

, oooooco 

. 0 . t * 
- 0. ' 

3.8000000 

.605874 

l .00000 

. 6000000 ' 

3 . 6000000 

l .21619 

1.00000 

3.4000000', 

r O'. 3- -7. :i - 

3 ♦ 400000O 

I . 86023 

1 .0000., 

3 . 2000000 

0 . 

3.2000000 

2.60063 

1 .OOOoO 

3 ♦ OOOOOOv 

0 . 

3 . 0000000 

3.60000 

1 . 000^ </ 

2.8000000 

0. 

2.8000000 

5.45391 

1.000 0<; 

2 . 6000000 

0. 

2.6000000 

13.5200 

1.00000 

2.5000000 

0. 

2.5000000 

0. 

1 .00000 


Choosing the gain to be 2.2 yields a closed loop transfer function as follows: 
Note that this matches the hand calculation. 


CLOSED-LOOP TRANSFER FUNCTION 


CLK= < Cl.NK/CLDK 2.200 


( 

( 


CLTF(S.) 

Cl-NPOLY < I ) 

2.200 )S ** 1 

5.500 ) 


NUMERATOR 

CL ZERO (I > 

( -2.500 ) + J< 0. 

CLNK*- 2.200 


) 


< 

( 

( 

( 

< 


CLDPOLY < I ) 

1.000 )S*t 4 < 
<6.000 )S** 3 < 
12.40 )S** 2 < 
13.20 )S** 1 < 
5 . 500 ) 


1 .034 
l . O'M 

) + Ji 0. 

-3.303 ) + J( 0. 

ci r>i<« 1.000 


CLTF(S) DENOMINATOR 

CLPOL E < I ) 
-.9039 ) + J< 

-.9089 > f j-; 

-.3790 


) 

) 

) 

) 


The unit step response follows as does the frequency response. 
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!’ 0 > 1 o i 

REGION DF RALCULAT JON: T» 0. 

A f- < T > = 0. 

TO V *» 
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I-' I 'v'0\ 5 1 ■; . 1 1,1 ,S I M 3 1 h t ) I,; JnXNI'U r / ‘JO i|l ‘IJ l ) v 


cons r ! 

PoIC --O.'JJ ' 

CONSTL = L .OP 
COiTi'1'2 - 13.73? 

HR AO •= 10.0? 

OELUCITY = 200. 0? 

CUiV:>f4 -- L.O? 

CQNST5 = 2.00? 

TAU1 - 0.40? 
r AU2 ~ 0.25? 

TAU3 = 0.30? 

TFINAL = 10.0? 

N ••= 5? 

PlllLIM •- 1.37? <* ENlER 
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IN RADIANS t ) 


WORD = PACKED ARRAY Cl.. 103 OF CHAR? 
RANGE — ARRAY Cl. *4J OF REAL? 


02200 ABSCISE A - ARRAY Cl.. 1001 3 OF REAL.? 

0 2300 ORDINATE - ARRAY Cl . . 1001 ? 1 . .23 OF REAL? 


02100 OAR 

02500 EPS , LEADX r LEADXD * DEL T , BE LTD * PHI * PHID » PHIDD » F'SI , PSID » T i 

02600 0 : ARRAY Cl.. 103 OF REAL? 

02700 I J ARRAY C1..10rl..23 OF REAL? 

02000 RANGE 1 J RANGE? 

02700 PLorrsx ; ordinate ? 



03000 

03100 

03200 

03300 

03400 

03500 


TIME ; ABSCISSA? 

I Y » P y M J INTEGER? 

I OPT t ICR i INTEGER? 
PRINTI J INTEGER? 

X » J» K J INTEGER? 
XCH» I CHAR J CHAR? 


; \> 


03600 

03700 

03300 

03900 

04000 

04100 

04200 


PLQ TIME J INTEGER? 
NLETTERS: INTEGER? 
TITLE t XWORD t YWQRD 
DELTAT J REAL? 

INC J INTEGER? 


WORD? 


. ■.) 

04300 

04400 

04500 

BEGIN <* MAIN PROGRAM * 

0 

04600 

TITLE := 'PSI OF A/C 


04700 

XWORD J= 'TIME 

i 

1 -v 

04800 

YWORD J= 'PSI 

' J 

04900 

T := o.o? 

\ 

05000 

LEADXD { = 0.0? 

• 

05100 

DELTD J= 0.0? 

r s) 

05200 

PHID := 0.0? 


05300 

PHIDD := 0,0? 

■ 0 

05400 

PLOTIME J= 1? 

05500 

psid :=o.o? 

;• O 

05600 

PRINTI := 0? 

05700 

LEADX 1= 0.0? 


05800 

CELT := 0,0? 

1 

05900 

06000 

phi := o.o? 

NLETTERS J= 10? 

0 

» 

06100 

xch := 'X'? 

[ ' 

06200 

06300 

IOPT .*=0? 
DELTAT :» 0.01? 

06400 

IY {=101? 

i. »• . _> 1 ' t 

06500 

P S* Is' i S 


REAL. ? 
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n , - i. » 

rOR K ; l TO 1 DO A 0 

RANGE ICKT 1= O.Ot C 

PS I 0 * 0 r 0 

WHILE T -> TFINAL DO 
BEGIN 

OK SET THE DER CVATI VE EGIJA ! J 
EPS : PStC - CuN3T4*P5I ••COr 

LEA.'ixn : (EPS-KAQNSTl- LEA O' 

BELT LEADX ! LLADXD# TAU I ? 
PH I! Ul.1 {-= <DELT*00NST2 -PHIDj ' 
PS ID J* (GRAV/ VELOCITY) *PIII. 
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: «) 

. M 'S ID 
AU2 i 


OK SHUNT LIMITS THE AIRCPAi i ROLL AH' X.E * ) 

IF (PIII.KPHID > 0.0) AND i^i.u'Hi) . i : L:_ LN ) fHLN PillD 0 : 0,0 


ASSIGN :JO 


/ -o k r» O 
\ ‘A J 

IC1 7 13 
ICi» 23 
I C2 f IT 
I C2 » 2T 
IC3f1T 
I C 3 » 2 T 
LZAtll 


I C5 f 1 T : = 


LEADX i 
LEADXD 
BELT J 
DELTD i 
PH I Li ? 
PH I BO t 

phi ; 

PHIDf 
PS I i 


PLQTPSICPLOTIMEf IT J= PS! 5 
I C5 f 2T :-= PSID i 

OK INTEGRATION *> 

FOR J := 1 TO N DO M 

ocjt:= icjfIT + icjfuj kdeltat; ^ 

OK REASSIGNMENT *) 

LEADX :« OCITf 
phid:= oc3t; 

PHI := 0C4T» 

psi := ocst; 

<* PRINT OUT THE PARAMETERS EVERY <BELTAT:K2S) INTEGRATIONS 


PRINTI J* PRINTI + 1? 

IF PRINTI = 1 THEN 
BEGIN 

WRITELNJ 

FOR X J= I TO N DO 

WRITELN<TfICXf1TfICXf2TfOCXT) : 

uriteln; 

PRINTI J= PRINTI - 25 i 

end; 

TIMECPLGTIMET J= Ti 
PLOTIME PL0TIME + 1 i 

t := t + deltat; 
end; 

ICHAR :■ 'X'f 

inc := i; 

WRITELN< ■'SUCCESS' ) » 


. r'r • 


X -V'T 
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TRACK I’NG,rEQUAT IONS 

jw-'i 'il. 


1. POSITION PARAMETERS: 



II HT || = R 


O (t) = ^ T (t) - iJj h 

-v 

/ 

/ 


(T/l#0£ •' v 


Pi. /■} ) , j 

\ 

\ O' 


f 

\ 


\ 1 
\ 

X 

' \ H rft / < fJ 


^ (Herr 


Let be the speed of the projectile fired on the 

pursuing host aircraft. The time taken to reach the target 

plane is T = V^/R. During this amount of time the target plane 

moves of an angle 6 equal to i^.T. That means that the 

projectile must be fired before ~ exactly when o = 6 . 

If we assume that a tracking system capable of giving the angular 

position of the target plane relatively to the host plane is 

available, the angular speed ij> can be computed as: 

(t+At) " ^ T (t) 

( t ) = - for every step of time At. 
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2. BASIC RELATIONS: 


I 
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In order to shoot the target as soon as possible, 
that is the next step of time, the next value of must 

ri 

be : 


(t + At) - ^ H (t) = <P H (t) X At 


A^ h + 6 = ip T ( t ) - ^ H (t) + $ T (t) x At 

v y / 

C7(t) ( Al^ T 



(t + At) 
^ R (t + At) 


t^ T (t + At) - 5 = ^ T (t) + iJ» T ( t ) x At - 6 
tjj H (t + At) - ^ H (t) 


At 
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We now need to generate a sequence of target positions, 
that is create the functions R and t|> T . Two ways were discussed: 

generating using the host plane parameters: 

the behavior of the target is conditioned by 
the behavior of the chaser. This is an approach 
very close to reality. 

programming functions and R in an appropriate 

way. We assume that the planes are similar in 

performance, particularly in speed, so that the 

range between them can only differ of R of a small 

value AR. t|; T must be a continuous function from 

TR + K 0, 2 IT ) . The end of fight is determined 

by T , that in real case can be interpreted 
max 

as a fuel shortage. 

The last option was chosen for our project, because 
it allows us to implement easily (change of function) a different 
target behavior. 
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CALCULATION ALGORITHM 


BEGIN. . 

INIT ip T (t) , R(t) , <P T (t) ; 

6 (t) = (t) x (V m /R(t); 

GET INPUT \p (t ) ; 

n 

/*Calculation of next ip: */ 


Do until TMAX: 

i^ H (t + At) = ip T (t) + | T (t) x At - 6 (t) ; (store ip^; ) 
(if * H (t + At) iP H max then * H = ifr Hmax ; & SE 


GET fi T (t + At) , R(t + At) ; 

. + At) = ^ T (t) . 

^ T (t + At) = ^ ' 

(if $ (t = At) ij/ T max then ip T = ^max;) 


6(t + At) = $(t + At) x (V m /R (t + At) ; 


NEXT: 

iKptt) ■*- ip T (t + At) 
4 > t ( t) <- * T (t + At) 
<5 (t) «- 6 (t + At) 


END; 


- '•> -.-Uto/go 4i 
• ,«\V* c ; 


ORIGINAL PAGE ES 
OF POOR QUALITY 


EXAMPLES OF TARGET MOTION GENERATION 


Tj%(t) = 2 it cos 2 tt t/t 

r T / max 

2 ir sin 2 tt t/t 

max 



(must be continued) 


R ( t ) = R (ip (t ) ) = a e (ellipse) 

a (1 + cos^) (cardioid) 

a = constant 


(must be continued) 

AR(t) = periodic or aleatory function giving small 
values (max - 10% . R) . 

(Must not be continued as it is a small value.) 
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LIMITATIONS 


In order to stay close to reality, we introduced 
the following limitations: The host plane turning rate 

cannot be greater than ( /20) radians/s, that is 9 degrees per 
second; the initial value of the target turning rate is set to 
0.01 radians/s, that is 0.57 degrees per second. 


CORRESPONDENCE ALGORITHM PROGRAM: 


Algorithm Program 


(t) 


OPSIT 

R 


RANGE 

$ T (t + 

At) 

NPSIT 

6 


DELTA 

At 


DELTAT 

*K (t + 

At) 

HPSID 

*H (t) 


PSI 

*H (t + 

At) 

EXPPSI 
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v *\ i. r/« - y f. 

In the computation of the tracking equations and 
anticipation parameters, converging results have been obtained 
with basic versions. 

The "fire" order will be displayed each timejijj 

Different test cases and air fight scenarios have 
been programmed and the results are satisfactory (fire orders 
displayed within 10 seconds — see listings attached) . 
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On the graph ,the positions of the host and the target 
airplanes are plotted . Each number on the points represents 
a time step . The + marks are the series of predicted tracking 


angles placed with the base of the next attaching the arrow 
of the previous tracking vector . This allows us to see 
that the series of optimal angles are drawing a curve looking 


like the target trajectory . This is a verification of the 
tracking angle generation procedure . The points 6, 7,8,9 


of the host plane trajectory show the optimal angle relative 
to the actual position of the plane . This is the real case. 
For each position of the aircraft the procedure gives 


the optimal angle the pilot has to command in order to shoot 
the target as soon as possible . As we can see on the graph , 
this command makes the plane point directly to the target . 


The firing condition is determined by a "window" of 
+ or - delta , delta being the lead angle depending on 
the range between the planes and the missile velocity . 
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NC f LF 7 

Cate 

lH-Uec-62 

it :2t:bu 

Monitor: 

use 

vOb 

FK735 

Rec 

ft 35 76 

f c r 

NETLFT 

Cate 

ltwJec-82 

l e :26:5b 

Monitor : 

use 

Job 

PHC T35 

Req 

<13576 

f o r 

NETLFT 

Cate 

let-Dec-82 
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u:2/:o2 Monitor: use 

•Start* 
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for NE UF7 

Cate 
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03576 
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FHCT35 Feo 

tt 35 76 

for NETLFT 

Cate 
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•Start* 

Job 
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0 35 76 

for NEUFT 
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CCCGOOOCCOOOOOCCCCOOCCCOOOOOOOOCCOCCCCCOCCOCCOOUCOCOCOOOOCCOCOOOCOOOOCOO 

rrfrnnnrrni i i i 1 i i 1 1 i i;u:93999v » it 



4^ 


ORlGlNAlj PAGE IS 
^QF POoS QUALITY 


o i rri 5'b i t u t. ‘<rz 4'u c < avci « ivt t rb^'OTTTrOTTiTTin*Tc*?TT?TrTTTTT7T , Tr^^^TBTT3TuT3 


CFI-CT^ Hecorclnc Irltlatec 
CL ink frcir BLUE N » TTY 303 


Set 18-Cec-E2 ^irOw'f'J 


TGFS-20 Comiranc processor 4(560) 
STY SIMCL.PAS 

PROGRAM SIM300F<INPUTfCUTPUT) i 


CONST 1 ’ 

FSIC =0.51 

consti =i.o; 

CCNS72 = 13.75 i 
C-RAV = io.o; 

, VELOCITY = 2C0.0! 
CONST 4 = 1.0 ; 
CONST f = 2. CO ; 
TAU1 = 0.40? 

TAL2 = 0.25! 


T A U 3 = 
Tf INAL 
N = 5 ! 
PHILIP 


0.50! 

= 10 . 0 ! 


= 1.57 


FI = 3.141552 
KISSIV = 1CCC 


<* ENTCR IN FAC1ANS O 


; VAR 


EPS *LEAOX.L£AO>Ot CELT ♦OELTC.FhI t FMC»FHCC»FiI,PSIC»T : REAL! 

C : ARRAY Cl.. 10] OF REAL! 

I .* ARRAY Cl.. 10,]. .23 OF REAL! 

FRINTI J INTEGER! 

>»J,K : INTEGCF! 

>CH , I CHAR J CHAR! 

R , CEL T AT : REAL! 

INC : INTEGER! 

CP3IT ,NPSI7 .RANGE , PS I TO , CELT A , E X PP S I ,H F S IC , E x h F HI , HP P I C ! REAL! 
PFSI C MX : FEAL ? 


BEGIN (* PAIN FRCGRAP *) 
T I = C.O! 

LEAOXC := C.O! 

CELTC != 0.0! 

FPTQ := C.O! 

fh ioc := o.o; 

FSIO !=0.0! 

FRINTI : = 0! 

LEAOX := 0.0! 

celt ;= o.o; 
f h i : r o.o ; 

CELT n := 0.10 ! 

RANGE != 2000.0! 


CPSIT != 0.5! 


FS ITC := 

C"LT/:= PS ITC*MSSLV/ RANGE! 

fs i := o.o ; 
hpsicpax := f i #2 o ; 

*H IC t ”"T~<=~TTTTmT _ DT' 

BEGIN 


K A IN LLOF 


l 



WRITELN? 

RZ= RANGE* ( 1*CCS (3PSIT > >/2 J 
WRITELN ( ‘RANGE = »»R)5 

writeln; 

£ F S := FSXC - CCNST4 * FSI - CCNS75*FS1C; 

LCADXO := (EPS^CCNSTl- LEACX)/TAL2» 

0 E 17 := LE A EX «LEAOXD*TAll ; 

P F IOO := (DEL7*CCNST2-FH IC >/7 AL3 J . 

PSIC := (GRAY/ VELOCITY) *PHi 

WRITELN (TTY , * ACTUAL HOST TLRMNG ANGLE : • » PSI ) ; 

W R ITELN 5 

(* SHUNT LIMTS THE AIRCRAFT RCLL ANGLE *> 

IF (PHl*FHIC > 0.0) ANO (ABS<FHI> >FHIL.1F) THEN PHIC : = Q.OJ 
(* CALCULATION OF THE EXPECTED FSI OF THE HOST PLANt. *> 
expp$i:= opsit«fs-tto*deltat; 

URITELN (TTY , ‘EXPECTED TURNING ANGLE : ‘ ) i 
(* SETS EXFFSI FCCULC 2«FI *) 

WHILE EXFFSI > 2 * FI DO EXFFSI 1 = EXFFSI - 2*FI ; 

WHILE EXFFSI < -<2«FI) CC EXFPSI *.= EXFFSI ♦ 2*FI ? 

URITELN (TTY , > *♦ EXFFSI) 5 

URITELN? - 

) If /eS(CFSIT - FSI) <= DELTA then 
/V WRITELN ( •<<..*** FIRE *** FIRE*** FIRt ***..5>')J 


oQ 


hfsic:= (Exffsi - psi > / celtat; 

URITELN ? 

UhlLE ABS(HFSIC) 5 FFSIOHAX CC 

HP SIC 1= (HFSICMAX * HFSIC / A8S(HFS1C>> 5 
WRITELN ( » ' ESTIHATED HOST TURNING RATE J • , HFSIL 
UR ITELNJ 


) ; 



(* 


GENEFATICN CF THE MOTION CF THE TARGtT FLANc. 




NPSITJr 2 « F I - ( 2»PI»CCS(2*FIMT*CELTAT)/7FlNAL> >? 

WRITELN ( * ANTICIPATED TARGET LINE OF SIGhT » * NFS1 T ) ; 

UR ITE LN J 


(* 


EXFPUSS HFS 3 D* WE LCCITY/GFAV; 
HFHIC:= (EXFFHI-FH > / CELTAT? 
FSITC:= (NFSIT-OPSIT) /celtat; 
celta:= ars(F3]tc«m:sslu/r > ; 

UR ITE LN (‘ANTICIPATED TOLERANCt 
UKITELN (* 

CPSIT := NPSIT ; 


(* EXPECTED VxLLE OF FH - ) 

(* COMMAND 10 BE GtNEHATE'D *> 

<* TLRMNG KATE OF THE TARGET * ) 
(- NEU ANGLt CF FIRING *> 

CcLTA ANGLE = »»CELTA)J 

> ; 

<* REASSIGNMENT «> 



■si 

IC 1,1] 


LEAD > i 

ic i»2 : 

* — 

LEAD >D J 

IC 2,1 3 

• z 

CELT ; 

ICS, 2 ] 

J z 

CELT c ; 

iCT,n 

J z 

FHID ; 

I C 3 f 2 1 

• z 

FHIDC ; 

I C A , 1 3 


f hi ; 

I C A , 2 3 

Z * 

fh c ; 



' •',* ' ■.V.'MfK* 
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) 



T C5.n := fsi ; 

w Cf,23 := fsid ; 

(* INTEGRATION *) 

FCF J := no f cc 

CCj:.*= I C C * 1 : 4 ILJ»2 3*CELTAI i 

(* REASSIGNMENT •) 

L S ACX := cc 1 3 ; 

Phic:= cc 3 ] ; 
ppi : = oca 3 ; 
psi := ccs 3; 

<* PRINT CUT TPE f ARA PETERS EVERY { CELT AT« 25 > IM t G PAT l CNS »> 
PP INTI := FFINTI 4 i; 

IF FRINTI = 1 THEN 
BEGIN 

WR I TEL N ? 

FOR X J: MO N DO ' ' 

UFITcLMT»ICX*13 tIC>t2 3*OCx 3) » 

WFITELN? 

WRITELN C ■ HCSl AIkCKAFT*) 5 

WRITELN ? 

URI T ELM 'EXFPSI = »»t>PPSIt» FPSIC = »*FFSID>5 

WRITELN? 

WFITELN ( *EXPPhI =- • »S >FPh I* » PFhlC = »»PFHO>? 

wpiteln; 

wkiteln; 

WRITELN < * TARGET * ) ? 

WRITELN? 

WFITELN ( *PSITD = • ,FS ITD ) ? 

WFITELNC ».\F$IT = • »NF S IT > 5 
hft ITELN ? 

WRITELN.? 

WRITELN ? 

PRINTi:= PRINT! - 25? 
lNC? 

t := t 4 celtat; 
en c ; 

WRITELN (• SUCCESS* ) ? 

ENC# 

a 

a 

a 

*■» 

SEX SIMCL.PAS 
LINK? Loacing 
CLNKxCT SIM3DC execution3 
INPUT : 

OUTPUT : 

CINFUTt end with *Z : 3 
"2 
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TIME STEF I 


RANGE = 


0.00 COCOE4 CC 

1 • 6775 £3c *03 


O.CU ( C00L40 C 


I 


ACTUAL PCST TURNING ANGLE 
EXFECTEC TURNING ANGLE : 


ORIGINAL 
OF POOR 


i 


PAGE 8S 
QUALITY 


l 



> 5.01CCC0E. .1 


ESTIMATED HOST TURNING RATE : 
ANTICIPATED TARGET LINE OF SIGHT 
ANTICIPATED TOLERANCE CELTA ANGLE 


1.57 07 96E-U 1 
1. 2396A8E- C2 
2. 5S69b4E-*-00 


C , n f' roOQE^C 0 
o . n u u o o o E • v •*. 
0.0C CC00E + C0 

o.oocoooc+oo 

O.OCCO0C£-»CO 


c. ocoooor^oo 

f. C 0 0 OOOE-O 1 

c . c • '-“ + 00 
o*cojoou r >oo 
c. coooooooc 


2,OGOCOOE*OU 
O.OOOOOOE+OU 
2. 2"C CCCE* 01 
C > u i. C. .j 0 E. ♦ 0 U 
C. GQCCOCE+OU 


2 • U C 0 LC OF - 01 
8.0C0CC0E-01 
2 . 2 C U l C 0 r -r o 0 

c.ococcof-*oo 

c . o c c c c c f ♦ o o 


HG3T AIRCRAFT 


EXFFSI 


5* 01 CO 0 OE-O 1 HfSIO 


1 .570 t S6E-01 


EXFFH I 


3.1«J59 iE^OO hFh IC 


3.141SS2L+Q1 


PSITD = 
NFS IT = 


TARGET 

-4.676C15E*C0 

1.233E48L-02 


TIKE STEF i 1.QQCC00E-C1 

RANGE = 1 • 9 S99 23E « 03 

ACTUAL HOST TURNING ANGLE : 0 . 0 0 0 C 0 OE* 0 C 

EXFEC7F0 TURNING ANGLE I 

> -<i. 7520 3CE-U 1 

<<..**• FIRE *** FIRE*** FIRE •**.♦0) 

ESTIMATED HOST 1URMNG RATE : -1 .57C79GE-U 1 

ANTICIPATED TARGET LINE OF SIGHT 4 .55A 4S1E- 02 

ANTICIPATED TOLERANCE CELTA ANGLc. - l.cS73«8t-01 


TIME STEF : 2.00CCC0E-CI 

RANGE = 1 • 9587 73E*03 


ACTLAL HCST TURNING ANGLE 
EXFECTEC TUNNING ANGLE : 


C,0uCC00f»0C 


•L 

O 


•3 


r i 



ORIGINAL: PAtSE.^ ls 
OFy PQQR v Q{jjALhTY3,'. 


— — > 8.G6F114E 2 

I ! ^ 

< < . . * V* FIRE *** FIRE*#* FIRE ***••>> 

ESTIMATED HOST TURNING RATE : , 1 .5 7C7 FGE- U 1 

ANTICIPATED TARGET LINE OF SIGhT 1 . 1 12 92 «E- C 1 

ANTICIPATED TOLERANCE CELTA ANGLE = 0. C E <=2 f 3E- U 1 


TIME STEF I 3.00 CC00E-C1 

RANGE = 1 • 9 938 J 3E ♦ 03 

AGTLAL HOST TURNING ANGLE : I.luOCOOE-03 

EXPECTED TURNING ANGIE : 

— > 1.73C399E-U 1 

<<••*** FIRE *** FIR£*«* FIRE ***••>> 

ESTIMATED HOST 1URMNG RATE : - 1..5 7 C 7 9GE* U 1 

ANTICIPATED TARGET LINE OF SIGhT 1. 973977E-CI 

; ANTICIPATED TOLERANCE CELTA ANGLE = A • 31 £bii 6E- 0 1 


TIME STEF I A.0QIC00E-C1 

RANGE = 1 • S £0 5 £ QE ♦ 03 

ACTLAL HOOT TURNING ANGLE I A.OlUOOOE-O 

EXPECTED TURNING ANGLE : 

— > 2.335031E-U1 

<<••*** FIRE *** FIRE*** FIRE ***••>> 

ESTIMATED HOST 1URMNG RATE I 1.57G796E-U1 

ANTICIPATED TARGET LINE OF SIGhT 3.075207E-C1 

ANTICIPATED TOLERANCE CELTA ANGLL = E. S 60 l*f OE- 0 1' 


TIKE STEF : D.jCCCCOr-n 

RANGE = 1.9S30£7E*03 

ACTLAL hCST TLRNING ANGLE I 


9.21 8CC0L-Q 


f W I <•# 
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i 


EXFE £ C TURNING ANGLE : 

> 4.1764:5E-Ul 

<<••*** FIRE *** FIRE*** FIRE ***..)> 

ESTIMATED HOST TURNING RATE J 1.57C7SGE-U1 

ANTICIPATED TARGET LINE OF SIGHT 4.412270E-C1 

ANTICIPATED TOLERANCE CEL7A ANGLE = 6 . £ <15 8 y 4E- 0 1 


TIME STEF : 6.00 CCOOE-C 1 

RANGE = 1.SC4228E+03 

AGTLAL HOST TLFNING ANGLE : 1.6y7718£-02 

EAFECTlC TURNING ANGLE : 

> 5.74F322E-U1 

<<•***• FIRE *** FIRE*** FIRE **♦«•>> ■ 

ESTIMATED HOST TURNING RATE t 1 .5 707 S6£-U 1 

ANTICIPATED TARGET LINE OF SIGHT 5.S7SES0E-Q1 

ANTICIPATED TCLERANCc DELTA ANGLE = £.2223uS£-01 


TIME STEF : 7.00CCQQE-C1 1 

RANGE = 1 • 8264 £9E 03 

ACTLAL HCST TURNING ANGLE : 2. 7.A 9279E-U 2 

EXPECTED TURNING ANGLE : 

— 7.S47509E-U1 

<<••*•* FIRE *** FI R £ * * * FIRE ***..5> 

ESTIMATED HOST TURNING RAT" : 1 • 57 C 7 S GE-u 1 

ANTICIPATED TARCt, L I N r OF SIGHT 7.7 71 675E-01 

ANTICIPATED TOLERANCE CELTA ANGLE = S . £ 11 1** 9E-0 1 


TIME STEP : 
RANGE = 
ACTLAL HCST 


6.0OCCCOE-C1 

1 .7123 6S£*03 

I 


A 


TURNING ANGLE : 


7* CF Alt-02 


. GRiGirf AL PAGE ;s 

• * ' l; ' 6F- P00R1QUALITY 



EXF. /EO TURNING ANGLE 


> 


S.563861E-U 1 


<<••*** FIRE *** FIRE*** FIRE ***«•>> 

ESTIMATED HOST TURNING RHE : 1 . 5 7 07 S6L- U 1 

ANTICIPATED TARGET LINE OF SIGHT S.7fcl J60E-01 

ANTICIPATED TGL £ R AN C £ CELTA ANGLE = ‘ 1. 1730 JSE* 00 



TIME STEF : 9.00CC COE-t 1 

RANC-E = 1.5S85£EE*03 

AGTLAL HOST TUFMNG'ANGLE I 4. 80 8C35E-0 2 

EXPECTED TURNING ANGLE : 

— — — — 1.17S04«E*U0 


<<••*** FIRE *** FIRE*** FIRE **-*.*>> 
ESTIMATED HGS T TURNING RATE : 
ANTICIPATED TARGET LINE OF SIGH 
ANTICIPATED TOLERANCE CELTA ANGLE = 


TIME STEF : - l.OOCCQOE+CO 

RANGE = 1.3G2376t*03 

ACTUAL HCST TURNING ANGLE ; 

EXPECTED TURNING ANGLE : 

<<••*** FIRE *** FIRE*** FIRE ****•>> 

estimated host turning rate : 
ANTICIPATED TARGET LINE OF SIGHT 
ANTICIPATED TOLERANCE CELT A ANGlL = 


1.5707S6E-U1 S 

t 

1.1S8SE1E+00 I 

1 ■ 42 35u It * 0 0 jj 

< 

i 

i 

f 

* i 

i 

5.835229L-02 
1.421846E*U0 

1.5707SEE-U 1 
l.441S07t>00 

1. 77t7faGE*00 



TIME STEP I l.lOCCOOE>OC 

RANGE - 1.1285 I3E*Q3 



L 


I 
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ACTIJ PCS T TURNING ANGLE : 6.F '423E-02 

V . U 

EXPECTED TURNING ANGLE : ' 

> 1.6638 I-3E-U G 

<<••»** FIRE *** F IRF*** FIRE ***..;> 

ESTIMATED HOST TURMNG RATE : 1.570796E-U1 

ANTICIPATED TARGET LINE OF SIGHT 1 . 702 5 2S£ - CO 

ANTICIPATED TCLZRANCc CELTA ANGLE = 2.313024E-CC 


TIME STEF I 1 .20 CO COE- CG 

RANGE = 6.6 62 4 13E- 02 . 

ACTLAL HOST TURNING ANGLE : 7.8«9617E-02 

EXPECTED TUkNING ANGLE : 

- — — — > 1 . S62972E-U C 

<<•.*** FIRE *** FIRE*** FIRE ***..5> 

ESTIMATED riCST TURMNG RATE : 1.57C7S6E-U1 

ANTICIPATED TARGET LINE OF SIGHT 1.9e2 046E-00 

U ANTICIPATED TOLERANCE CELTA ANGLE = 3 . 2 1 4 6 4 CE ♦ 0 0 

TIME STEP I 1 .30 CC COE- CC 

RANGE = G.GC2423E-02 

ACTUAL hCST TURNING ANGLE l E. 9166116-02 

EXFECTEC TURNING ANGLE : 

> 2 .2611 E6E-U 0 

<<..♦** FIRE •** FIRE*** FIRE «**»•>> 

ESTIMATED HOST TURNING RATE J 1.S70796E-U1 

ANTICIPATED TARGET LINE OF SIGHT 2.273121E-C0 

ANTICIPATED TOLERANCE CELTA ANGLE = 4.9327J.6E-00 


7 1 TIME STEP : 1 .4 0 CO 00 E- C 0 

. L . 1 

RANGE = 


3. 5019 COE- 02 
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i 

ACUA^CST TURNING ANGLE : , S.Ss«/C05E-0J 

EXPECTED TURNING ANGLE : . 

— > 2.574214E+U0 

<<••*** FIRE * * * FIRE*** FIRE ***••>> 

ESTIMATED HOST TUR MNG RATI : 1.5707S6E-U1 

ANTICIPATED TARGET LIRE OF SIGHT 2.5SQQ20E*CO 

ANTICIPATED TOLERANCt CELTA ANGLE = t . S 062W 5c *-0 0 


TIME STEF : 1.50 £ COOL* CC 

RANGE = 1.4£29£'6E*02 

;AC T LAL HOST TURNING ANGLE : 1 • 0 y 7 J 20E- C 1 

EXPECTED TURNING ANGLE J 

> , 2.9Qiyo?E*U 0 

<<••*** FIRE *** FIRE*** FIRE ***••>> 

ESTIMATED HOST TURNING RATE : 1.57Q7=6c~Ul 

ANTICIPATED TARGET LINE OF SIGHT 2 . 9 16 4 fi Ac ♦ C 0 

ANTICIPATED TOLERANCE CELTA ANGLE = 2. 2014U QE *01 


TIKE STEP J '1 .6 00000 E* CO * 

RANGE = 2.5230C5E*01 

ACTUAL HOST TURNING ANGLE 5 l.ly5e35E-Ul 

EXPECTED TURNING ANGLE I 

— — — — --> 3.24294FE*U0 

<<••«** FIRE -** FIRE*** FIRE ***..>> 

ESTIMATED HOST TUNING RaTF : 1.57G7F6E-U1 

ANTICIPATED TARGET LINE OF SIGHT 3. „ .,o 2 36c* CO 

ANTICIPATED TOLERANCE CELTA ANGLE = 1. 2 4661 4t *02 


TIME STEF 


1 .7 0 £ 0 0 OE* 0 C 
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R ^ .N G t 6.5643 '^q'^T 

AC TLAL HCST 7LFN ING ANGLE { 1.3U2555E-01 

EXPECT EC TURNING ANGLE : 

— > 3.5$5987E*uG 

<4..*** FIRE *** FIRE*** FIRE »**••>> 

ESTI.’UiuO tW:T ’"«MNG RATE : 1.57C7S6E-U1 

ANTICIPATED TARGET LINE OF SXGhT j , cO / 5 33E* CO 

ANTICIPATED TOLER AN CL CELIA ANGLE = f . I b ? 7 u 5L ♦ 0 2 


TIME STEF : 1 .6 0 £ 0 G OE* C C * 

RANGE = I • C £78 C2L* 02 

ACTLAL HOST TURNING ANGLE : 1. 4u 52 /8L-0 1 

EXPECTED TURNING ANGLE : 

> 3. 5556 30E*U0 

<<••*** FIRE **• FIRE«** FIRE ***.*)> 

ESTIMATED HGST TURNING RATE : 1 .5 70756E-U 1 

ANTICIPATED TARGET LINE OF SIGHT 3. 970 1 cSE* 00 

ANTICIPATED TOLERANCE CELTA ANGLE = 3. 3 525 o 3c ♦ C 1 


TIME STEP : 1,90 COOOE+CO 

RANGE = 3.2408 £ 3c ♦ 02 

ACTUAL HOST TURNING ANGLE : 1.5U7SS8E-01 

EXPECTED TURNING ANGLE : 

---> 4.3324 43E-MJ 0 

<<••*** FIRE FIRE*** FIRE »*♦.»>> 

ESTIMATED HOST TURNING RATE : 1 . b 7 C 7 8 6 L- U 1 

ANTICIPATED TARGET LINE OF SIGHT 4. 3 41 572E* 00 

ANTICIPATED TOLERANCE CELTA ANGLE = 1 . 14 593 4c* G 1 


TIME STEP 


2.00 C000E+C0 


• <* 
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RANGE G.3762 26t>02. ^ 

ACTUAL HOST TURNING ANGLE : ■ ; . 1.610717E-01 

EXPECTED TURNING A N G L £ : 

> 4.712955E+U0 

<<••*** FIRE * * » FIRE*** FIRE ***••>> 

ESTIMATED HOST TURNING RATE : 1.57C7S6E-U1 

ANTICIPATED TARGET LINE OF SIGHT A . 7 2 0 6 1 6E* 0 0 

ANTICIPATED TOLERANCE DELTA ANGLE = S • 5 A 46 / 7t* 0 0 


TIME STEF I 2.10CGOOE*CO 

RANGE = 1 • Q Cfl 2 « 9t •» 03 

ACTLAL HOST TURNING ANGLE : W 7 i 3 A 26c.- 0 1 

EXFECTED TURNING ANGLE : 

.... — — > 5.Q9«:6G<i£*UC 

<<...*. FIRE *** FIRE*** FIRE ***..>> 


ESTIMATED HCST 1UPMNG RATE : 
ANTICIPATED TARGET LINE OF SIGHT 
ANTICIPATED TOLE fiANC i CELTA ANGLE = 


1.5707S6E-U1 
5. 1 C5€ 2 IC + CC 
2. E206*s 1E*00 


TIME STEF : 2.20 CQOOE* CC 

RANGE = 1 « 2 £3 3 £9E* 03 

ACTLAL HCST TURNING ANGLE : 1.8I61SPE-01 

EXFECTEC TURNING ANGLE : 

> 5 • 4 S 10 RAE + U 0 


ESTIMATED HOST TURNING RATE J 1.570756E-U1 

ANTICIPATED TARGET LINE OF SIGhT 5.A556S0E*0Q 

ANTICIPATED TCLtRANCc CELTA ANGLE = 2.tlPln6t*0C 


» * 


TIME STEF 


2.30CO OOC- CO 
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i 

R*NC:.~- 1 . 7C56 2 2c. *1)3 , " 

ACTUAL HOST TURNING ANGL't* 1 -; *' 1* 91 c B 75E- 0 1 

EXPECTED IUkNING ANGLE : 

> 5.3A554SL+U0 


ESTIMATED HOST TURNING RATE : 
ANTICIPATED TARGET LINE OF SIGHT 
ANTICIPATED TOLERANCE CELTA ANGLE 


1 .57C7SGE-U 1 
5 • t So fc 57E ♦ 00 
2. 3 03 9t> 3E*0 0 


TIME STEF I 2.40 (CODE* CO---- 

RANGE = 1.52317e£*u3 

«CTIAL HOST TURNING ANGLE : 

EXPECTED TURNING ANGIE : 



estimated hcst turning rate : 

ANTICIPATED TARGET LINE OF SIGHT 
ANTICIPATED TGLERANCE CCLTA ANGL t = 


2 . U 2 j.j 9 SL -0 1 

6.2E1625E+UQ 

1.S7C7S6E-U1 
G.2 6 3131E + 00 
2. 05141 be ♦ 0 0 


TIME STEP : 2.50 COOOZ* CC 

RANGC = 2.0000 C 0 £ ♦ 0 3 

ACTUAL HCST TURNING ANGLE : 2 . 1 * 4 3 1 4E - U 1 

EXPECTED TURNING ANGLE : 

. > 3.545217E-U1 


ESTIMATED riOST T U R M N « " ' : 

ANTICIPATED TARGET LINE OF SIGHT 
ANTICIPATED TGLERANCE CELTA ANGLE = 


1 ,M7fi7S6E-U 1 
6.G / 7 70 GE ♦ CO 
1. S7262l£«-00 


2.50COO 0E*00 
2.5C COOOS + GO 
2.5C COO CZ ♦ 0 0 
2 . 5 C C C 0 0 E ♦ C G 
2.5 C COO CE ♦0 0 


1.07 4 179 E- 01 
6.(7169^-02 

c. ouonooooo 

2.G54383Z+0C 
2. 12431 4E-C1 


-1.027524E-01 
C.COCCOCE-*OU 
1 • 4 t 9 70 6£ *0 U 
C.COCCOCE+OU 
1.027194E- 01 


5.754266E-02 
fc.e 1 1£5 4? -02 
1 .43 c 97CbF-01 
2.0542E3E + 00 
2.227C22E-01 
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We have successfully demonstrated the ability 
of the On Board Simulator and Tracking Procedure in high 
level PASCAL code . The implementation details have yet 
to be completely worked out . 

The On Board Simulator was tested using a st,ep 
input and studying its response . The aircraft turning 
time constant is 4.4 seconds and damping ratio is .65 , 
which provided satisfactory results . The final simulator 
will fly at a velocity of 1000 feet per second and have 
a maximum bank angle of 82.8 degrees for an 8g turn ( without 
aerodynamic limitation ) . 

The angle tracking is done relative to an inertial 
reference angle stored in the host aircraft's reference 
system . Quick encounters are simulated by inputting a 
programmed target maneuver into the host aircraft's field 
of view . 

A more accurate tracking method could be developed 
by following the host and target aircraft separately in an 
inertal grid and computing relative information from the 
inertial systems . This would allow longer time engagements . 

The On Board Simulator would not be easily upgraded 
from its present 3 degrees-of-freedom . This would require 
additional graphics on the screen display and would require 
much additional work on generating relative dynamics . 
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I. System Description 
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A. Purpose 

The objective of the project is to design a micro- 
processor based Aircraft Spin Warning System which periodically 
samples the assymetric thrust and yav; rate of an airplane and 
then issues voice synthesized warnings and/or suggestions to 
the pilot of how to response to the situation. 

The system is to meet the requirements set forth in 
the June-August 1982 status report of the system study for the 
Application of Microcomputers to Research Flight Test Techniques, 
as summarized in table 1 of the paper (included in this 
report in Figure 1 ). 

Such a system is expected to aid the pilot in recovery 
from spins and high speed departures which occur during flight 
tests of aircraft at flight boundary limits. 

B. Type of Microcomputer to be Used 

Our S'«S design is based around' the Motorola MC68000 
16-bit microprocessor, implemented on a M68KVM01AZ monoboard 
microcomputer. Since the SVJS is a real time application where 
computing speed is critical, a very fast microprocessor must 
be selected. The advanced contemporary design, 8-MHZ clock 
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Figure 1* System Block Diagram 
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rate, and 16 -bit word width all contribute to fast processing 
times and make the MC68000 a good choice for this application. 

Convenient I/O and bus structures provide additional 
advantages that simplify and smooth system design and integration 
since we have access to peripheral devices with compatable I/O 
buses. 

C. Selection of Speed Synthesizer Device 

The SWS project is a developmental project and so 
requires versatile design features to allow for the many design 
changes that always occur during the development of a system. 

This is the main reason behind the selection of the SC01 from 
the various voice synthesis products on the market. The 
phoneme based SC01 allows us to speak any word desired, and 
easily accomodates any changes in the future. LPC based speech 
methods provide better voice quality, but have restricted and 
unversaitile vocabularies and are much more expensive. The SC01, 
on the other hand, has been on the market for a couple of years 
and has well documented dictionaries and application articles 
making it a clear choice for a short term, low budget project. 

D. Voice Synthesis Design Description 

The block diagram of the voice synthesis design is 
shown in Figure 1 . The one byte of discrete data, seven 
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channels of ADC data, and PTM timer are controled by the input 
interface. The two byte output data (phoneme and pitch) is 
transferred to SC01 voice synthesizer via the output interface 
as well as to the hexadecimal LED Display for test of the data. 

The software moniter aircraft parameters from the ADC 
and discrete input register, analyzes these parameters and 
then, if necessary, sends voice synthesizer to speak out the 
required message. 

The functional details of the modules and their design 
are described in the next two chapter. 
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II. Software Development 


A. General Description 

As indicated in Table 1, there are a total of lit 
sentences to be spoken in various cases. These sentences 
are 'UNLOAD', ' LEFT-ENG-MIL ' , ' RIGHT-ENG-MIL' , ' LEFT-ENG- 
IDLE ' , ' RIGHT- ENG- IDLE ' , ' BOTH- ENG-IDLE ' , ' STICK-FULL-LE^T ' , 

" STICK-FULL-RIGHT ' , ' STICK-FULL-FWD ' , ' CENTER- CONTROLS ' , 
'RECOVERY-COMPLETE', ' STICK-HALF-AFT-AND-HOLD ' , 'LEFT-ENG- 
OFF', ' RIGHT- ENG-OFF' . 'In some conditions, depending on 
input values received from ADC channels, up to 3 sentences 
need to be spoken. Each sentence may contains 8 to 26 
phonemes . For each phonemes , the SC01 takes an average of 
100 ms to complete its speech synthesis process. Therefore, 
in some cases, the time consumed to speak three sentences 
may take more than 2 seconds. Two seconds is a large 
amount time for computer busy waiting for phoneme output. 

Note that another requirement in Table 1 is to keep sensing 
each ADC channel to determine if the situation has changed 
or not. The sensing rate should be much less than 2 seconds. 
It is obvious that there are two requirements to be per- 
formed: (l) the software has to keep sensing ADC channels 
while the SC01 is speaking some sentences. (2) a real time 
clock is necessary to count the actual time that the 
temperature exceeded the temperature tolerance limit or for 
how longer the YAW RATE has recovered to its normal 
condition, etc. 

Therefore, a program 'VOICE' and a subprogram ' 
SMALLMOUTH' are implemented to solve the above requirements. 
The program 'VOICE' is used to dump sentences in an output 
buffer, the subprogram SMALLMOUTH will send all phoneme and 
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pitchcod.es from the output buffer to the SC01 via 2 8-bit 
ports for voice synthesis. These two programs interleave the 
CPU by timer and parallel port interrupt . The data structure 
for the sentence table and buffers are shown in Figure 
Each buffer is a record which contains a bufferfull flag, 
along with 50 phoneme codes and 50 pitch codes . The program 
VOICE contains a procedure SET- SENTENCE- TABLE which stores 
the phoneme and pitch symbols of all sentences. Then VOICE 
sets up the timer clock, senses ADC channels when timer 
interrupts arrive, and then according to the ADC channel 
value, decides which sentences are to be spoken and dumps 
them into output-buffer. The subprogram will send phoneme and 
pitch codes from the outputbuffer to the SC01 and wait for 
the phoneme to finish its sound generation. After d phoneme 
code has completed its sound generation process in the SC01, 
the SC01 will send back a parallel port interrupt request to 
execute the subprogram SMALLMOUTH again and get next phoneme 
and pitch code. 


The program VOICE is clussfied under the following 
headings: (l) Initialization procedures. These include the 
procedures to set up the sentence table, to set up the real- 
time interrupt clock, to set initial conditions and allow 
interrupts. (2) a MAIN body. This includes a procedure 
NEW-CASE which resets software message flags, two external 
procedures GET-DISCRETE and GET-VALUE to get the discrete 
values from parallel port: #2 Jand get YAW RATE, ALTITUDE, IAS, 
A0A,AY , AIRSPEED and TEMPERATURE from ADC channels. After 
obtaining those ADC vuluds, procedure MAIN which contains 
many CASE blocks, will select the sentences to be spoken, 
a procedure BIGMOUTH-SPEAK will copy those sentences from the 
sentence tuble into buffers in translation -queue. 
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A procedure SPEAK is responsible for dumping the sentence 
buffers in the translation-queue into the output buffer. An 
external function OUTPUT-BUFFER-EMPTY will check the buffer-full 
flag of both output buffers. If any of the output buffers are 
empty, the procedure LOAD -OUTPUT-BUFFER will load a sentence 
buffer from translation-queue into the output buffer. The way 
to load a sentence buffer to the output buffer is as following: 
First find which of the sentence buffers in the translation- 
queue needs to be dumped , then convert the phoneme and pitch 
in the sentence buffer to ASCII code, then an external proce- 
dure DUMP will dump the phoneme and pitch codes into the 
output buffer. After the whole sentence is loaded into the 
output buffer, an external procedure. -MARK-OUTPUT-BUFFER-FULL 
will set the buffer-full flag of the output buffer and also 
wakes up the subprogram SMALLMOUTH. The subprogram SMALLMOUTH 
will send the phoneme and pitch codes from the output buffer 
one by one to the SC01 for voice synthesis. 

After the procedure SPEAK dumps all the sentences into 
the output buffer, the program VOICE completes its job in this 
time frame and then the program goes into an external procedure 
HALT. The procedure HALT has a loop body. The CPU will keep 
looping in this loop body until the next timer interrupt comes. 

B. Abstract Data Structure 

An abstract data type is given as follows: 

Abstract Data Type 

Sentence-buffer = Record 

Buffer-full flag 

Phoneme symbol array [ 1 . . 50 ] 

Pitch symbol array [ 1 . . 50 ] 

End 


Translation cjueue = Array [ 0 ..2 ] of sentence buffer. 
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Output-Buffer = Record 

Buffer- full flag 

ASCII phoneme code [ 1 . . 50 ] 

Pitch code [ 1 . . 50 ] 

End 


Abstract Program Structure 
Declare 

VOICE ( 0 ) > OUTPUTBUFFER 

SMALLMOUTH ( OUTPUTBUFFER ) — Voice generation 


Begin 


When timer-interrupt comes 
do VOICE 

When parallel port interrupt comes 
do SMALLMOUTH 


End 


Program VOICE 
Begin 

INITIALIZATION ; 

Repeat 

MAIN (* Get ADC value, decide sentences *) 
SPEAK (* Dump sentence to output buffer *) 
HALT (* Wait for next timer interrupt *) 
Until FOREVER (* or turned off *) 


End 
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Subprogram SMALLMOUTH 
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Begin 


Declare Waked up 
if output buffer is empty 
then return interrupt 
else 

send out pitch code 
send out phoneme code 
increase pointer 

endif 


End 


C. Implementation 

The first level of the abstract data structure imple- 
mentation contains four procedures: INITIALIZATION, MAIN, SPEAK, 
and HALT. The process SMALLMOUTH is also implemented. The 
SMALLMOUTH contains a set of instruction to send out phoneme and 
pitch codes to the parallel port for the SC01 to speak. This 
subprogram is called up by a parallel port interrupt or a trap 
from VOICE. SMALLMOUTH send out the next phoneme along with 
pitch data, then returns control to VOICE. Under normal 
conditions, no warning signal need be generated. VOICE finds 
nothing to speak, then resets the runningflag and halts. 

SMALLMOUTH recognizes the situation from the runningflag and halts 
also. The whole system is in an inactive situation, only the 
procedure VOICE will be called by the timer every half second. 

If the situation does not require a voice message, the system 
becomes inactive again. 

A more detailed description of the different levels of 
implementation is shown in Figure 2 and the attached program 
listing. 


INITIALIZATION MAIN SPEAK 
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D. Detail Explanation of Each Procedure 

1. Program VOICE. 

The program VOICE will dump the sentence need to be 

spoken into out put buffer as follow: 

Step 1. a) initialize the timer clock as 0.5 second by calling 
external procedure SET -REALTIME- INTERRUPT. 

b) set a look up sentence table by calling internal 
procedure SET-SENTENCE-TABLE. 

c) set initial condition for variables by calling 
internal procedure SET-INITIAL-CONDITION. 

d) allow interrupts by calling external procedure 
ALLOW- INTERRUPTS . 

Step 2. a) initialize a new case condition when a new timer 
interrupt comes by calling internal procedure 
NEW-CASE. 

b) get the value of discrete by calling external 
procedure GET-DISCRETE. 

c) get value of YAW RATE, ALTITUDE, IAS, AIRSPEED, 

AOA, TEMP. 

find the sentence and pass this sentence as a 
parameter to internal procedure BIGMOUTH- SPEAK 
and load the sentence into a sentence buffer in 
the translation- queue by calling this procedure. 

Step 3. dump all the sentence buffers into output buffer. 

Step 4. wait for next timer interrupt. 

Step 5. go to step 2. 

Step 6. end. 


2. Procedure INITIALIZATION 

This procedure sets the initial conditions such as 
sentence table, PTM clock rate, etc. This procedure will be 
executed only when the system is turned on. 
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Step 1. set realtime interrupt. 

Step 2. set sentence table. 

Step 3. set initial condition 
Step 4. allow interrupts. 

3. Procedure SET-REALTIME- INTERRUPT 

This routine sets the cascading timeout interrupt of 
the MC6840 PTM, connect 03 to C2 and get realtime interrupt 
from timer //2. It also sets the realtime interrupt clock 
from clock §1 of PTM. The timer will wake the procedure MAIN 
every half second, get ADC value, and decide what sentences are 
required to speak. 

4. Procedure SET-SENTENCE- TABLE 

This procedure sets up a sentence table where all phoneme 
and pitch symbol are stored in l4 SOUND-BUFFER as a look up 
table. 

5 . Procedure SET-INITIAL-CONDITION 

This routine sets the initial condition for variables 
TEMP-COUNT, COUNT-30, COUNT-40, and flags such as WARNINGFLAG, 
FOUND, etc. 

Step 1. reset all software synchronization flags 

Step 2. call external procedure MEMORYINITIALIZATION to 
initialize the variables and flags in assembly 
part of the programs. 

6. Procedure ALLOW- INTERRUPTS 

When power is turned on, all interrupts are disabled. 

This assembly routine store a $0480 in the control register 
of M68K000 versa system which will enable the timer and 
parrallel port interrupt. Also, this routine store autovector 
addresses for timer interrupt, parrallel port interrupt, and 
trap #6. 
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T. Procedure MEMORYINITIALIZATION 

This routine is called by the procedure of SET-INITIAL- 
CONDITION to reserve block space in memory for output buffers. 
It also clear buffer-full flags of output buffers , set ARBIT 
point to the buffer #0, declare that the subprogram SMALLMOUTH 
in not running. 


Step 1. a) ARBIT: reserve one byte for ARBIT, ARBIT 

point to the output buffer is dumping. 

b) SLEEP: reserve one byte for SLEEP, SLEEP 

flag to indicate the SMALLMOUTH is not 
active . 


c) BFOFUL : reserve one byte for BFOFUL, BFOFUL 

flag to indicate that buffer 0 is full. 

d) BF1FUL: reserve one byte for BF1FUL, BF1FUL 

flag to indicate that buffer 1 is full. 

e) BUFO: reserve 50 word to store phoneme and 

pitch codes. 

f) BUF1 : reserve 50 word to store phoneme and 

pitch codes. 

g) QUEPTR : reserve one byte for QUEPTR, QJJEPTR 

point to the offset in output- 
buffer is speaking. 


Step 2. 

ARBIT 

//0 

: set ARBIT to point to output- 
buffer $0 . 

Step 3. 

SPK 

ttO : 

set speak point to output buffer HO 

Step 4. 

QUEPTR 

#0: 

point to the first word in output 
buffer . 

Step 5* 

SLEEP 

#1: 

declare that the SMALLMOUTH is not 
active. 

Step 6. 

BFOFUL 

//0: 

reset buffer-full flag of output 
buffer if 0 

Step 7- 

BF1FUL 

m-. 

reset buffer-full flag of output 
buffer ft 1. 


Step 8. 


End. 


X rt» 
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8. Procedure MAIN 

The procedure MAIN will get input conditions such as 
YAWRATE,, AOA, TEMPERATURE, IAS, AIRSPEED, AY, ALTITUDE, from 
ADC channels and decide which sentences need to be spoken. 

Then copy the sentences from the sentence table into buffers in 
translation- queue. Another procedure SPEAK will move the 
sentences in these buffers into an outbuffer which is located 
at some relocatable address in memory. Then the subprogram 
SMALLMOUTH will move each phoneme and pitch codes to the SC01. 


Step 1. declare new case. 
Step 2. get DISCRETE. 

Step 3. get value of 

a. IAS. 

b. AIRSPEED. 

c . AOA . 

d. AY. 

e. TEMPERATURE. 

f. YAWRATE. 

g. ALTITUDE. 


Step’’4. determine YAWRATE is positive or negative. 

Step 5- if TEMPERATURE > 121 5°C for 3 second , BIGMOUTH 
speak 'LEFT OR RIGHT ENG OFF'. 

Step 6. if AOA < -10° and AY > .lg, BIGMOUTH speak ' 
STICK HALF AFT AND HOLD ' . 

then if A0A>0 for 5 second, BIGMOUTH speak 1 
RECOVERY COMPLETE'. 


Step 7* if warning flag = true 

then begin if YAWRATE <1 40°/sec for 2 seconds, 
BIGMOUTH speak 'CENTER CONTROLS', 
if YAWRATE < 30°/ sec for 2 seconds , 
and AIRSPEED > 120 KTS, BIGMOUTH 

speak 'RECOVERY COMPLETE'. 


end 


: . . 5 . r s J.- i-X Jw 
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Step 8. if any condition in Table 1 occured, then 
BIGMOUTH speak sentences according Table 
and set warning flag. 

9- Procedure NEW CASE 

This procedure initialize a new case condition when 
a new timer interrupt comes. 

Step 1. reset all buffer-full flag in three sentence 
table of translation queue. 

Step 2. call procedure MEMORY-RESET to initialize; a 

new case condition in assembly part, also clear 
interrupt flags. 

10. Procedure MEMORY-RESET 

This procedure will be merged in the procedure of 
NEW - CASE to reset the initial condition and forget the old 
sentence in the buffer when a new case happens. It first 
clears the interrupt request flag and then allows all timer 
to operation mode by writting #0 in control register 1 ( 

CR20 = 0 ). Next, it clears PTM status register by clearing 
interrupt request flag. It continue by reading timer #2 counter, 
then by testing ARBIT which is pointer to one of the output 
buffers, then it clears the output buffer which ARBIT is 
pointing to. It returns after reseting port interrupt flag. 

11. Function GET-DISCRETE: integer 

This function get the value of discrete from parallel 
port and return the value as an integer. 

12. Function GET -VALUE 

The function GET-VALUE converts the variable name to 
channel number, gets the ADC value from the forward function 
ADC, then returns the channel value to the specified variable. 


) 
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Step 1. convert variable name to channel number. 
Step 2. get ADC value from the according channel. 


13. Function ADC 


The function ADC will return the ADC value from the 
specified channel. 


Step 1. 

Step 2. 
Step 3. 
Step 4. 
Step 5« 
Step 6. 
Step 7- 
Step 8. 

Step 9* 

Step 10. 


set address of base, control register, MSB 
and LSB. 

load control register to sleet channel number. 

initial conversion. 

wait until ready. 

move MSB to data register #0. 

shift data register #0 left 8 bit. 

move LSB to data register HO . 

arithmetic shift data register HO right 4 bit to 
get 2 ' s implement value . 

store data register HO to functional value 
return address. 

return. 


14. Procedure BIGMOUTH-SPEAK 

This procedure gets a sentence name from the procedure 
MAIN, then copies the sentence from sentence table into a 
sentence buffer in translation queue. 

Step 1. copy sentence from sentence table to sentence 
buffer which is pointed by NEXTEMPTYBUFFER . 

Step 2. set NEXTEMPTYBUFFER point to the next sentence 
buffer. 


15- Procedure COPY 

This procedure copies each phoneme and pitch symbol from 
one sentence is sentence table to a sentence buffer in 
translation- queue. After all phoneme and pitch symbol are 
copied, the buffer- full flag of the sentence buffer is set. 
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16. Procedure SPEAK 

This procedure sends out sentences from translation 
queue to output buffer whenever the output buffer is empty. 

After all sentences have been sent out, the queue is empty. 

Step 1. find if any one output buffer is empty. 

Step 2. then load output buffer 

Step 3. if translation queue is not empty, go Step 1. 

17. Function OUTPUT-BUFFER-EMPTY: boolean 

This function responds with a boolean value to indicate 
if there is an empty output buffer. If it finds an empty 
output buffer, it returns a function value of true. This function 
starts by finding if ARBIT points to BUF1 or BUFO. The next 
stepis to find out if the buffer-full flag of either BUFO or 

A 

BUFl is set. If it is not set then it returns the value of 
true. 


18. Procedure LOAD-OUTPUT-BUFFER 


This procedure dumps the sentence from the translation- 
queue into output buffer, each phoneme and pitch symbol will 
be converted into ASCII code and then dumped into the output 
buffer . 


Step 1. 
Step 2. 

Step 3- 
Step L. 
Step 5- 


find next sentence-buffer ready to speak. 

a) convert phoneme code. 

b) convert pitch code. 

dump to output buffer, 
mark output buffer full. 

reset the buffer-full flag of the sentence buffer. 


19. Function CONVERT-PHONEME-CODE 

This function convert phoneme symbols into ASCII code 
before dumping them into the output buffer. 
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20. Function CONVERT-PITCH-CODE 

This function converts pitch codes from decimal to 
octal before the pitch is dumped into output buffer. The 
first digit in decimal form will send to SC01 for pitch control. 
The next two digits will be used as frequency control for the 
SC01 clock input. 

21. Procedure DUMP 

This procedure dumps pitch and phoneme codes into the 
empty output buffer at the offset location INT. The pitch 
code is dumped in the lower byte, the phoneme is dumped in 
the higher byte of the word. The procedure starts by getting 
the offset INT, pitch code and phoneme code. Then it test 
ARBIT to see which buffer it is pointing to. If it is pointing 
to BUF1 then load buffer #1 and if it is pointing to buffer it 0 
then load buffer #0. 

Step 1. a) get offset INT. 

b) get pitch code 

c) get phoneme code. 

Step 2. test ARBIT. 

Step 3. if ARBIT = 1, goto step 5. 

Step 4. load base address of BUFO to register A2, goto 
step 6. 

Step i>. load base address of BUF1 to register A2. 

Step 6. a) move pitch code to buffer, 

b) move phoneme code to buffer. 

Step 7- end. 

22. Procedure MARK-OUTPUT-BUFFER-FULL 

This procedure is called at end of the procedure of LOAD- 
OUTPUT-BUFFER when the whole sentence in translation buffer have 
been dumped into output buffer. The buffer-full flag of the 
output buffer which is pointed by the ARBIT will be set. Then, 
this procedure will check the flag SLEEP to determine that if 
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the subprogram SMALLMOUTH is in active or not. If the SLEEP 
flag is set , go trap #6 to wake the subprogram which will send 
phoneme and pitch codes to parallel port #1. 

23. Function QUEUE- EMPTY: boolean 

This function checks the three sentence buffers in the 
translation-queue. If all buffers in the queue are empty, it 
returns a true variable. 

2k. Procedure HALT 

After execution is completed, this procedure is called 
in the main program that cause the main program to become 
idle. It first resets the interrupt request flag to wait for 
next timer interrupt. Then, it keep checking to see if there 
is an interrupt from the timer and then it exits after timer 
interrupt and it also sets the flag after the timer interrupt 
arrives . 

Step 1. reset flag to wait for timer interrupt. 

Step 2. test 1 timer has interrupted ?' 

Step 3- if timer has not interrupted go to step 2. 

Step 4. a) exit after timer interrupt, 
b) set flag. 

Step 5* return from interrupt. 

Step 6. end. 

25 • Subprogram SMALLMOUTH 

This subprogram is initialized by parallel port interrupt 
or trap #6. The parallel port #1 is used to output phoneme and 
pitch codes. When a phoneme is completely spoken by the SC01, 
the SC01 returns an interrupt request via the lower byte of #1 
parallel port. It will initialize this subprogram to send out 
the next phoneme and pitch. Another chance to execute this 
subprogram is after a sentence is dumped into an output buffer. 
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the procedure MARK-OUTPUT-BUFFER-FULL will check the flag 
SLEEP and generate a #6 trap which autovector address is the 
starting address of SMALLMOUTH. After a sentence is completely 
sent out from output buffer to SC01, the flag SPK is changed to 
indicate that another output buffer is the next to be send out. 
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F. Code Generation 

The software including all assembly and pascal routine, 
was developed and debugged on the TQP20 system. Next, all 
the assembly and pascal routines were transferred to IBM U3U1 
system by operator of the USC Engineering Computer Lab. The 
purpose of transferring our software from TQP20 to IBM U3U1 
system was to sue the MC68000 support software package on IBM 
system to get the object codes for all assembly and pascal 
routine and link them together. 

a. MC68000 Support Software Package. 

This package can be used as a powerful software develop- 
ment tool for MC68000 based system. The package utilities enable 
the production of relocatable machine code for the Motorola 
MC68000 microprocessor. The source code can be written in pascal 
or the assembly language of the microprocessor. A linkage 
editor can create execution modules in which library functions 
can be selectively included. The package provides libraries 
for floating point operations as well as for runtime routines 
for three popular MC68000 based configurations including versa- 
raodule, the one used for this project. 

The support software is designed to operate in a host 
system. In this case, the host is an IBM ix3Ul called VIRGIL. 


25 
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VIRGIL is part of the USC Engineering Computer Lab. It has 
8 mega bytes of main memory and it runs under CMS (conver- 
sational machine system). 

b. Pascal Cross Compiler 

Processes pascal programs and produces relocatable code 
that can be subsequently linked with other modules. The 
compiler works in two stages. The first one, phase 1, checks 
the syntax and semantics of the source code and produces and 
interned iate pcode file. The second stage, phase 2., processes 
this intermediate code and produces a relocatable and position 
independent object code module. 

c. M68000 Cross Macro Assembler 

It processes 68000 assembly language files and produces 
object code files that can be subsequently linked with other 
files. 


After using the linker to produce the complete and 
final object code package, the code was down loaded to the 
lab's computerra computer (an Intel 8085 bs?ed CP/M microcomputer) 
Finally, the code is dumped into the 68000 for execution. 
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III. Hardware Development 

A. MC68000 System: 

/ 

The MC68000 nonoboard system with the PTM has the 
following capability: (1) Vector Interrupt Handler. (2) 
Programmable Timer module which enables the real time 
interrupt and adjustable pitch control for SC01 Voice 
Synthesizer. (3) The processor instruction set provides 
software interlocks for processes to interleave the CPU. 

Trap instruction and TAS (test and set) are useful for mul- 
tiprocess communication that support modem structured 
programming techniques. 

B. Input Interface 

law rate , Altitude, IAS, AY, Airspeed, AOA and Temperature 
are transfered via the M68MML5A 16 channel high-level A/D module. 
The 68000 processor also gets the discrete signals fromthe #2 
parallel port without handshaking. 

C. Output Interface 

Phoneme and pitch codes (two bytes) are transfered from 
the 68000 system to the SC01 Voice Synthesizer utilizing the 
l6bit capability of the parallel port. Design and layout of 
both the discrete fetch and Voice Synthesizer was completed as 
the attached schematic figure. The driver LS2ljL receives the 
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phoneme code from lower byte of parallel port #1. These phoneme 
lines will be pulled up by a U.7K resistor array, then sent to 
the CMOS SC01 device and to the displays. A strobe signal from 
pin P1CA2 of the parallel port interface will initiate the voice 
conversion. After the phoneme has completed, the SC01 will 
request the next phoneme code via PlCAl vector interrupt request. 
MC68000 system also provide a 2Mhz clock from FTM #3 output line 
to the 7ii97 Binary Rate Multiplier for pitch c ontrol. A 8-bit 
pitch control code will be provided by the upper byte of parallel 
port #1, which will be latched by LS37ii D-type transparent 
latches. Two MSB bits are sent to the SC01 Cor direct pitch 
control. The LSB 6 bits sent to the Binary Rate Multiplier will 
be able to adjust the SC01 input clock rate from 30K hz up to 2M 
hz which will provide a minor adjustment of pitch control. An 
amplifier LM386 will provide audio amplification to drive a 
speaker. A 10K ohm potentiometer provides volume control. 


i 
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Figure U. Board Layout 
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Figure 5. External Hardware Block Diagram 
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IV. Results / Conclusions 


A. Test Setup 

For test purposes, a modification of the routine to 
obtain data from the ADC board was made. Instead, switches on 
the external hardware board were used to simulate certain 
combinations of aircraft parameters. The switch pattern was 
input through the discrete data input port, and the system 
responds to this input. 

To initiate the system testing, a small assembly language 
test routine was used to test out the interface between the 
68COO and' the external hardware board. This routine sent the 
word "unload'' continuously to the SC01 display devices, and 
provide valuable in trouble shooting the interface. 

B. Results / Problems Encountered 

Successful operation of the system was never achieved. 

The last of a series of problems encountered was the inability 
of the linkage editor to provide the start address of the loaded 
object code. 

V/e were able to get the small test routine to get the system 
to speak "unload" for a short time, but then the SCC1 failed for 
no apparent reason. 

Work yet to be completed includes determining the start 
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address of the run-time package, system software checkout, 
voice pitch control adjustments, and system integration testing 
with sensors providing the ADC inputs. 
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1 SC-01 SPEECH SYNTHESIZER 

A Division of Federal Screw Works 


■- . . 

500 Stephenson Highway 


DATA SHEET 

Troy, Michigan 48084 

> . 



Votrax® CMOS Phoneme Speech Synthesizer 


GENERAL DESCRIPTION 

The SC 01 Speech Synthesizer is a completely self-contained 
solid state device. This single chip phonetically synthesizes 
continuous speech, of unlimited vocabulary, from low data 
rate inputs. Figure t. 

Speech is synthesized by combining phonemes (the building 
blocks of speech) in the appropriate sequence. The SC-01 
Speech Synthesizer contains 64 different phonemes which are 
accessed by a 6-bit code. It is the proper sequential 
combination of these phoneme codes that creates continuous 
speech. 

The SC-01 Speech Synthesizer is cost-effective, consumes 
minimal power and enables in-house product development 
without vendor dependency. Signals from the SC-01 are 
applied to an audio output device to amplify and distribute 
the synthesized speech. See Figure 2. 



Figure 1. Votrax® SC-01 Speech Synthesizer 
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• Single CMOS chip 

• 70 bits per second 

• 22 pm package 

• 9 ma. current drain 

• Wide voltage supply range 

• Latched 5V. compatible inputs 

• Digital pitch level inputs 

• Automatic inflection 

• On-chip master clock circuit 

• Optional external master clock 

• Variety of voice effects 

• Sound effects 

• Customer product security 


The design of the equipment specified herein is proprietary. 
Rights for the reproduction and distribution of the data contained 
herein are granted except for the manufacture and reproduction 
of the subiect equipment. 


v* Vp 



Figure 2. SC-01 Flow Diagram 


Votrax® reserves the right to alter its product line at any time, or change specifications or design without notice and without 
obligation. 
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PHYSlC’Al! DESCRIPTION 5 

The SC-01 Speech Synthesi/ei isa22 pm Large Scale Integrated 
Circuit which contains all the circuitry necessary to generate 
phonetically synthesized speech. The SC-01 is fabricated using 
CMOS technology, wmch offers high input impedance and low 
power drain 

ELECTRICAL DESCRIPTION 

The SC-01 Speech Synthesizer is a program-compatible with 
existing Votrax® phoneme synthesizers. It requires 70 bits ol 
data per second for continuous speech production. The 6-bit 
phoneme codes are 5 volt logic compatible and are latched for 
data bus applications A phoneme-construction algorithm and 
filters, within the chip, create the synthesized audio output. 
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PHONEME DESCRIPTION 

Table 1 lists the 64 phonemes produced by the SC-01. Each 
phoneme code is accompanied by its symbol, average duration 
time, and an example. The underlined segments of the 
example word demonstrate the phoneme use, i e., sound to be 
pronounced. 

Table 2 subdivides the 64 phoneme symbols into seven 
categories. Each category represents a different production 
feature. The first six categories are characterized by voiced, 
fricative (expired voice), and nasal sounds. The seventh 
category is characterized by phonemes with no sound output 

PHONEME PROGRAMMING 

Manual Operations: Votrax® maintains a library of 

phonetically programmed words. Reference to this library and 
programming manuals will aid in word synthesis. 

Automatic Operations: Votrax® can supply a micro-computer 
system for automatic conversion of English text into phoneme 
sequences. This system is particularly useful tor m-house 
vocabulary development and product secunty Contact 
Votrax® for further information. 
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Figure 3 SC-01 Footprint and Outline Dimensions 
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Table! Phoneme Chart 
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Phoneme 

Code 

Phoneme 

Symbol 

Duration 

(ms) 

Example 

Word 

00 

EH3 

59 

jacket 

01 

EH2 

71 

enlist 

02 

EH1 

121 

heavy 

03 

PA0 

47 

no sound 

04 

DT 

47 

butter 

05 

A2 

71 

made 

06 

A1 

103 

made 

07 

ZH 

90 

azure 

08 

AH2 

71 

honest 

09 

13 

55 

inhibit 

0A 

12 

80 

inhibit 

08 

11 

121 

inhibit 

0C 

M 

103 

mat 

00 

N 

80 

sun 

0E 

B 

71 

bag 

0F 

V 

71 

van 

10 

CH* 

71 

chip 

1 1 

SH 

121 

shop 

12 

Z 

71 

ZOO 

13 

AW1 

146 

lawful 

14 

NG 

121 

thing 

15 

AH1 

146 

father 

16 

001 

103 

looking 

17 

00. 

185 

book 

18 

L , 

103 

land 

19 

K 

80 

trick 

1 A 

J* 

47 

judge 

18 

H 

71 

hello 

1C 

G 

71 

get 

ID 

F 

103 

fast 

IE 

D 

55 

paid 

IF 

S 

90 

pass 


/T/ must precede /CH/ to produce CH sound. 
/D/ must precede /Jt to produce J sound. 


Phoneme 

Phoneme 

Duration 

Example 

Codej 

Symbol 

(ms) 

Word 

20 j 

A 

185 

day 

21 \ 

AY 

65 

day 

22 j 

23 .? 

Y1 

80 

yard 

UH3 

47 

mission 

24 5 

AH 

250 

mop 

25 | 

P 

103 

past 

26 } 

• 0 

185 

cold 

27 i 

1 

185 

pin 

28 ’ 

U 

185 

move 

29 ; 

Y 

103 

any 

2A 

T 

71 

tap 

2B 

R 

90 

red 

2C 

E 

185 

meet 

2D 

.W 

80 r .... 

win 

2E 

£E~ 

J8C.7; 

dad 

2F 

AE1- ' ‘ 

•> 103 1 

after 

30 

AW2 

90 

salty 

31 

UH2 

71 

about 

32 

UH1 

103 

uncle 

33 

UH 

185 

cup 

34 

02 

80 

for 

35 

01 

121 

aboard 

36 

IU 

59 

you 

37 

U1 

90 

you 

38 

THV 

80 

the 

39 

TH 

71 

thin 

3A 

ER 

146 

bird 

38 

EH 

185 

get 

3C 

El 

121 

be 

3D 

AW 

250 

call 

3E 

PA1 

185 

no sound 

3F 

STOP 

47 

no sound 




Table 2. Phoneme Categories According to Production Features 


Voiced 




'Voiced' 

Fricat. 

‘Voiced* 

Stop 

Fricative 

Stop 

Fricative 

Nasal 

No 

Sound 

E 

EH 

AE ‘ 

UH 

001 

Z 

B 

T 

S 

M 

PA0 

El 

EH1 

AE1 

UH1 

R 

ZH 

D 

DT 

SH 

N 

PA1 

Y 

EH2 

AH 

UH2 

ER 

J 

G 

K 

CH 

NG 

STOP 

Y1 

EH3 

AH1 

UH3 

L 

V 


P 

TH 



1 

A 

AH2 

0 

IU 

THV 



F 



11 

A1 

AW 

01 

U 




H 



12 

A2 

AW1 

02 

U1 







13 

AY 

AW2 

00 

w 
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SIGNAL DESCRIPTION (See Figures 4 and 5) 


NOTE 
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Phoneme 6-Bit Selection Code (P9-P5): Data input is to six 
pins. Latching is controlled by the strobe (STB) signal. 

Strobe (STB): Latching occurs on rising edge of strobe signal. 

Inflection Level Setting (1 1, 12) Instantaneously sets pitch level 
of voiced phonemes 

Acknowledge/Request (A/R): Acknowledges receipt of 

phoneme data (signal goes from high to low one master clock 
cycle following active edge of STB signal). Also indicates 
timing out of old phoneme concurrent with request for new 
phoneme data (signal goes Irom low to high) 

NOTE 

If external phoneme timing is desired, phoneme 
requests can be ignored. However, best speech is 
realized with internal timing. 

Master Clock Resistor-Capacitor (MCRC): This input 

determines the internal master clock frequency. Select R-C 
values for 720 kHz to achieve standard phoneme timing. 
Connect this input to MCX when using internal clock, ground 
when using external clock. 


Varying clock frequency varies voice and sound 
effects. As clock frequency decreases, audio 
frequency decreases and phoneme timing 
lengthens. Figures 6 and 7 illustrate manual and 
DAC (Digital to Analog Converter) voice variation 
schematics, respectively. 

Master Clock External (MCX): Allows control by an external 
dock signal. 

NOTE 

Ground MCRC during MCX operation. 

Audio Output (AO): Supplies analog signal to audio output 

> : ^ ffu,' 

ew '“‘ YP !A UQ fiP0 £ ! V) 

Audio Feedback (AF): Used with Class A or Class B transistor 
audio amplifiers for added stability. 

Class B (CB): Current source for Class B transistor audio 
amplifier. 


Table 3. Timing Specifications 


CHARACTERISTIC 

SYMBOL 

MIN 

TYP 

MAX 

UNIT 

Input Setup Time (P, to STB) 

T S 

450 



NS 

Input Hold Time (P ( to STB) 


» 



NS 

Rise Time of STB Edge (.8V to 4V) 

t rs 



109 

NS 

A/R Width (A/R Connected to STB) + 

t arw 

1 

1.3 

2 

MS 

STB Width 

T sw 

299 



NS 

STB Low* 

Tsl 

• 



NS 

Propagation Delay (STB toA/R after Tarw* 

t dar 



500 

NS 

A/R Rise Time (Capacitive load = 36pf) 

t rar 



109 

NS 

A/R Fall Time (Capacitive load = 3®pf ) 

t far 



109 

NS 

Time from A/R Request to STB Service) 

t ars 

9 

< 

500 

M* 

Time of Phoneme Duration * 

t ph 

47 

197 

259 

MS 


t Dependent on Master Clock frequency 720kHz 

’ Strobe must remain low (72x Master Clock Period) before rising edge 
Votrax® reserves the right to alter its product line at any time, or change specifications or design without notice and without 
obligation 
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j A/R CONNECTED TO STROBE 


QSIGINAL PAGE IS Figure 4. Timing Diagram 
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NOTE 



Figure 6. Variable Voice by Potentiometer Control 

Vp 



Data at address 0 must be a pause phoneme code. j| 

Multiple Message, Fixed Block Size: See Figure 9. Message 
address block is loaded into the counter. The message is then 
clocked out of the ROM by the A/R signal. 

NOTE 

Message Block = 2 n maximum. 

Multiple Message, Variable Block Size: See Figure 10. The 
microprocessor loads phonemes into a data bus. The A/R 
signal generates an interrupt request for each new phoneme. 

CONNECTING THE AUDIO OUTPUT DEVICE 

Audio Output: The AO signal has a max imum peak to peak volt 
age swing of .26 times Vp. depending upon the phoneme se- 
lected. and the AO signal is D.C biased 


Figure 7. Variable Voice by DAC Current Injection 

Class A Amplifier: See Figure 1 1 For a single transistor 
amplifier, the selection of R. C, or R s values depends upon the 
TYPICAL APPLICATIONS value of Vp and the desired audio level. 

General: The SC-01 Speech Synthesizer is easily designed into 
systems ranging in complexity from ROM/counters to 
microprocessor controllers. 

Single Message System: See Figure 8 When the counter is re- 
leased (START isTRUE), the message is clocked out of the ROM 
by the A/R signal. The system must be stopped when DONE is 
TRUE Note When using A/R tied to STB, connect a .01 uf 
capacitor to TP3 to insure power up reset of SC-01. 


Vp 


o 

COUNTER ROM SC-01 ^ R 



Figure 8. Single Message System 


Votrax^ reserves the right to alter its product lint at any time, or change specifications or design without notice and without 
obligation 
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Figure 9. Multiple Message, Fixed Block Sure 


Class B Amplifier: See Figure 12. A current source (CB) is 
required for this push-pull amplifier. 

NOTE 

Minimum power is consumed when speech is 
inactive. When Vp = <-12.0 volts and R s = 40 
ohms, the bias current drain is approximately 3.5 
milliamps. 

Controlling Audio Output Power: See Figure 13. A resistor or 
potentiometer from the speaker to ground can be used to 
control the audio output power. 



Figure 10. Multiple Message, Variable Block Sue 



Figure 1 1 . Class A Amplifier 


Figure 12. Class B Amplifier* 


Figure 13. Controlling Audio Output Power 


For Class B Amplifier: ( P) x (R s min.) - 81.6 x (Vp) where {3 is beta or current gain of transistor. The AO line is protected by an 
internal series current limiting resistor of 90 ohms maximum. If more current is required of the SC-01, then the above formula 
indicates distortion will occur. 


Votrax® reserves the right to alter its product line at any time, or change specifications or design without notice and without 
obligation. 


- :.■( 
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Table 4. Analog Output Specifications 



CHARACTERISTIC 

MIN 

MAX 

UNIT 

Output Voltage (AH Phoneme) 

.18 x Vp 

.26 x Vp 

Vpp 

Output Bias Current * * (.6V<C8<Vp) 

3.5 

7.3 

mA 


ELECTRICAL CHARACTERISTICS: T„ = 9 to 70°C, Vp = 7 to 14 V E 


' Variable 

“With CB. AF, AO connected t or Class 8 audio amplifier (see APPL ICA TION NO TES) 
* Frequency of Master Clock ^.1.25 / RC 
Note 1P1, TP2 must be left open f oi normal operation. 



CHARACTERISTIC 

MIN 

TYP 

MAX 

UNIT 

Digital Input Impedance 

1 meg 



Ohm 

Input Capacitance <Pj , STB) 



3 

P< 

Input Capacitance (11, 12, MCX) 



8 

P< 

Digital Input Logic "0" (except 11, 12. MCX) 

V G - 0.5 


< 

Cl 

* 

© 

00 

< 

o 

o 

Digital Input Logic "0" (MCX) 



V G * 1 0 

< 

o 

o 

Digital Input Logic "0" (11, 12) 



2 x Vp 

Vdc 

Digital input Logic "1" (except 11,12, MCX) 

V G + 4 0 


Vp • 0 5 

Vdc 

Digital Input Logic "1” (11, 12) 

.8 x Vp 



V DC 

Digital Input Logic “1" (MCX) 

4.6 



V DC 

Digital Output Logic "0" (1 sink = 0.8rnA) 



V g +0.5 

v DC 

Digital Output Logic "1" (1 source = 0.5mA) 

Vp— 0.5 



V DC 

Power Supply Current Vp - 9V 


9.1 


mA 

Vp = 9V" 


11 

18 

mA 

Vp = 14V" 


18 

27 

mA 

'Master Clock Frequency 


720K 


Hz 

MCX Input Duty Cycle 

60.40 


40:60 

% 

Master Clock Resistor Value (MCRC)”’ 

6.5k 



Ohm 

Master Clock Capacitor Value (MCRC)*** 



300 

pf 


Votrax® reserves the right to alter its product line at any time, or change specifications or design without notice and without 
obligation 
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Table 5. Absolute Maximum Ratings 
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OF POOR QUALITY 


A8SOLUTE MAXIMUM RATINGS * 


RATING 

SYMBOL i 

VALUE 

UNIT 

Power Supply Voltage 

Vp 

20 

Voc 

Power Dissipation at 25° C 

P 0M i 

I 

650 

rnW 

Derating Above 25° C 

l 

1 

5 

mW C 

Operating Ambient Temperature 

T ° ! 

0 to 70 

X 

Storage Temperature 

t stg ; 

--55 to 1 25 

’c 

Input Voltage 


-0.5 to Vp+0 5 

U 

O 

> 

DC Current Max. Above Vpr0.5V 

'iNM 

1 0 

ma 

Lead Temperature (soldering 10 sec.) 


300 

’c 


* Operation above these limits could damage the device. 

i 

NORMAL OPERATING CONDITIONS: 7v ^ Vp .$ I4v, Q° C £ T 0 $ ?Q°C J 




Votrax® reserves the right to alter its product line at any time, or change specifications or design without notice and without 
obligation. 
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JZ. 

PI PftO 

33 

P 1 P A 1 

31 

PI PftZ 

35 

P 1 P A 3 

33 

PI PRH 

3 1 

P 1 P R 5 

31 

PIPRId 

Z1 

P 1 P A 1 

35 


PICftZ 43 


CLOCKS* i 
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ITY IZONElLTR 


REVISIONS 


DESCRIPTION 



ft LDECI 


B RDEG 

C 

D 

VCC 

VLED 

BIRMK ■ 

STROBE 
G N D , 

noi-R • 


DATE APPROVED 
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OISPLM 


LQEC ay 

mc M 


ftUDIQ iW 

oo 


VCC 

VLED 

BLANK 

/STRQSE 

6 ND 

noi-K 


DISP.LAV 















+ IZVDC 

- -- GND 
























®UT t 
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0UT - 


PICA I 


PZPAO 
PZ P A I 
PZPftZ. 
PZPA3 
PZPAi 
PZPAS 


B 


GO 



t 


DWG NO 







FILE: TEST V5 


ASM 


A 


i 


VM/SP CONVERSATIONAL monitor system 


************* ****$******$$**$*****$******$#***************$*******>>:** 
* 

* 

** PROGRAM TEST 1 
* 

* 

* THE PROGRAM TEST1 OUTPUT A SENTENCE '’UNLOAD” TO SC01 

* WITHOUT USING INTERRUPT AND TIMMER. 



SECTION 

PHNPOT 

EQU 

PITPOT 

EQU 

PHNCTR 

EQU 

PITCTR 

EQU 

LATCH 

EOU 

START 

MOVE.B 
MOVE. B 
MOVE. B 
MOVE. B 


CLR.L 
CLR.L 
CLR. L 
MOVE A. T. 
MOVEA. I, 
RR A . S 

LOOP 

MOVE 

COUNT 

SUBQ 
RNE. S 


SUBQ 

BNE.S 

SPEAK 

MOVE. B 
MOVE. B 
MOVE.B 
MOVE- B 
MOVE. B 


CMPI. B 
BEQ. S 
ADDQ 
BRA 

EXIT 

ORA 

TABLE1 

DC 


DC 


DC 

TABLE? 

DC 

TABLE 3 

DC 

DC. L 


9 


SF70021 PORT #1 LONER BYTE 

TP70020 PORT #1 UPPER BYTE ^ 

SF70025 PORT #1 LOWER BYTE CONTROL REGISTER 

$F70024 PORT #1 UPPER BYTE CONTROL REGISTER 

TPF /- 


TABLE2 .PITPOT SEND OUT PITCH CODE 

W, PITCTR 

SLATCH, PITCTR SEND OUT LATCH SIGNAL 

#0 , PITuTR 


DO 
D 1 
D? 

STABLE 1 , A2 
STABLE!, A4 
SPEAK 
#570. D2 
# 1 .02 
COUNT 
# 1,01 
LOOP 

0 fA2, DO) .PHNPOT 

#0 . PHNCTR 

#L ATCH .PHNCTR 

#0, PHNCTR 

0 (A4,D0) , D 1 

#*?' ,0 (A2,D0) 

EXIT 

# 1 , DO 

LOOP 

START 

' 2MX57 • 

S 1 E 


$5B 

103 

71 


SET BASE ADDRESS TABI.E1 
SET BASE ADDRESS OF TABLE2 

SET INSIDE LOOP COUNT 
TNERT LOOP BODY 
KEEP LOOPING 
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PILE: TESTV5 ASM A VM/SP CONVERSATIONAL MONITOR SYSTEM 


DC • L 

103 

DC 

121 

DC. B 

90 

DC 

55 

DC. B 

127 

END 



APPEHDIX C 


AIR DATA SYSTEM 

SOFTWARE SOURCE CODE 


FILE: EQU2 


PASCAL 


A 1 


1 

I 


VM/SP CONVERSATIONAL MONITOR SYSTEM 


• PROGRAM MCH ; 

FUNCTION MI (PTI , PSI : REAL) : REAL ; 

CONST A,B,C,D,E,F,G; 

VAR X: REAL ; 

BEGIN 

X; - pTI/PSI : 

IF X <= 1-893 THEN 

MI := SQRT(5. 0) *SQRT (( (PTI/PSI) ** (2/7) ) -1) 

ELSE 

MI := ( (A*X-B*X-C*X*X-D*(X**3) -E* (X**4)-F* (X**5) -G*(X**6) -G* (X**9) ) /X) ; 
END: 

BEGIN 

END. 


FILE: EQU 3 PASCAL A1 VM/SP CONVERSATIONAL MONITOR SYSTEM 


subPROGRAM STATPR ; 

FUNCTION FPSINF (PTI ,MI:rEAL) : REAL; 

VAR A : RE A L ; 
begin 

A :=1 . 0/ (7 . 0* MI*MI) ; 

IF MI <= 1.0 THEN 

FPSINF := PTI/exp (3. 5*in ( 1 +0 . 2* MI* MI ) ) 

ELSE 

FPSINF := PTI*A*exp (2. 5* In (1 -A) ) /0. 1839371 ; 
END. 


FILE: EQU5 PASCAL A1 VM/SP CONVERSATIONAL MONITOR SYSTEM 

subPROGRAM OUCC : 

FUNCTION FQCC (PTI, PSI N F : R EAL) : R EA L; 

BEGIN 

FQCC:= PTI-PSINF; 

END. 
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FILL): EQU6 


PASCAL 


VM/SP CONVERSATIONAL MONITOR SYSTEM 


subPROGRAM KNOTS: 

FUNCTION FKEAS ( Si NF , PSI NF : REA L) ; REAL ; 
BEGI N 

FKEAS := MINF*66 1 . 48*SQRT (PSINF/2 1 1 6 . 22) ; 
END. 


FILE: EQ07 


PASCAL 


VM/SP CONVERSATIONAL MONITOR SYSTEM 


subPROGRAM CALSPEED; 

FUNCTION FKCAS (QCC: S E AL) : REAL ; 

BEGIN 

FKCAS := 1479. 1 *SQRT (exp (0.28571*ln (1+ QCC/ 21 16. 22) )-l) ; 
END. 


FILE: EQU8 


PASCAL 


VM/SP CONVERSATIONAL MONITOR SYSTEM 


subPROGRAM GEO: 

FUNCTION FHP (P$INF, PSI : REAL) : REAL ; 
VAR A, R: REAL ; 

BEGIN 


IF~R S > N o{ 2 2 3 36 I^EIEN FHP : = 1 4 54 42* (e xp (0 . 1902 632 365*ln (r) ) -1) 
ELSE IF R > 0.0540328 THEN FHP : = 164 21 9- 2080 5*LN (PSI) 


ELSE BEGIN 

A : = exp(0. 01463563358*1:1 (0.0540328/R) ) : 

IF R > 0.00856663 THEN FH P : =7 10 79 4* A* A- 6 45 1 77 
ELSE BEGIN 


END 

END; 


A : =exp (0.04097 977 *ln (0. 0 0 85666 3/R) ) ; 
FHP : = 8166 0.7* A* A- 16 29 28. 8; 


I 


FILE: ADSTIM 


ASH 


VM/SP CONVERSATIONAL MONITOR SYS 


* * *** * **♦ * * * * 4c * * * 4c * 4c * * * * * * * * ♦ * * * * * * * * * * * * * * * * * * ♦ * * * * * * * * * * * * * * * * * * * * * * 


ADSTIH ASH 
WHITTEN BY 
PURPOSE 


FADI KUBDAHI AND CABPO SOSA 

THESE ROUITNES ABE USED TO MEASURE THE EXECUTION 
OF EQUATIONS IMPLEMENTED IN ADS. THE TIMMING IS 
IN THE FOLLOWING MANNER: 

SETIHE : 1* INITIALIZE, LOAD AND START TIMER *) 
EQU#; (* INVOKE OPERATION TO BE TIMED *) 

COUNT := READTIM ; (* STOP TIBER, READ COUNTER 


* 

* 

TIME* 

DONE* 

* 

♦ 

* 

*) * 


WUULU « — unailiiu | V X A JL ULlU m ilijn U lliUU T l T 

******************* 44*4*44*44*444**444*»*44444*444l<4*44*444444**»***44 

*******1f^^*^if*********************** ********************************** 


ROUTINE 

PURPOSE 

FORMAT 


METHOD 


S£TIH E 

: INITIALIZE AND START TIMER 
: USED AS A PASCAL PROCEDURE 
PROCEDURE SETIME : 

: TIMER MODE IS CONTINOUS WITH 
TIMER COUNTER IS LOADED WITH 


#3 OF THE MONOBOARD’S 
DECALRED AS: 


PTM 


AN INPUT OF 
ITS MAXIMUM 


2 MHZ. 
VALUE 


_ (FFFF) 

** *** *********** *** ********* ** ******* * *************** *** * * *** «****« 

MSB BUFFER REGISTER (WRITE ADDRESS) 

LSB BUFFER REGISTER (READ ADDRESS) 
WRITE TIMER # 3 LATCHES 
WRITE TIMER # 2 LATCHES 
BEAD TIMER # 3 COUNTER 
WRITE CONTROL REGISTER # 

WRITE CONTROL REGISTER # 

WRITE CONTROL REGISTER # 

CONTINUOUS MODE CODE 
SINGLE SHOT MODE CODE 


MSBBUF 

EQU 

IF70009 

LSBBUF 

EQU 

SF7000F 

T3LTCH 

EQU 

$F7000F 

T2LTCH 

EQU 

SF7000B 

T3CONT 

EQU 

IF7000D 

CNTRL3 

EQU 

$F70001 

CNTRL2 

EQU 

SF70003 

CNTRL1 

EQU 

SF70001 

T3MODE 

EQU 

$80 

T2MODE 

EQU 

$A 0 


XDEF 

SETIME 


SECTION 

9 

SETIME 

EQU 

* 


3 (CR20=0) 
1 (CR20= 1) 


MOVE. L 
MOVE. B 
MOVE. B 
MOV E. B 
MOVE. B 
MOV E. B 
JMP 


i^i^DE. CNTRL3 
#1 ,CNTRL2 
#0 ,CNTRL 1 
t$FF , MSBBUF 
#$FF,T3LTCH 
(A4) 


********************************************************************** 

* ROUTINE READTIM * 

* PURPOSE : READ THE CONTENTS OF THE COUNTER OF TIMER #3 * 

* FORMAT : USED AS A PASCAL FUNCTION DECLARED AS: * 

* FUNCTION READTIME : INTEGER : * 

* OUTPUT : READTIM : CONTENTS OF THE COUNTER. * 


4 ********* ***** * ********* **** * **************************************** 


LI 


EQU 


READTIM 


XDEF 

SECTION 

EQU 

MOVE. L 
LINK 
MOVE. L 


READTIM 

9 

il‘J5Lr (,7 > 

A6, 12 (A5) 
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FILE: ADSTIM 


ASM 


VM/SP CONVERSATIONAL MONITOR SYS 


CLB.L 
MOVE. B 
LSI 

MOVE.B 

MOVE.! 

UNLK 

MOVE. L 

RTS 

END 


D 1 

T3CONT, D1 
#8, D 1 
LSBBUF.D1 
D1, 12 (A6) 

A6 

(A7) +, 12 (A5) 


I 


APPENDIX D 

USER'S DOCUMENTATION 
FOR SYSTEM SUPPORT SOFTUARE TOOLS 
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Communication between COMPUTER , 1 < CEE-560L) and USC-5CL CEngr 


Computer Lab) 


I . ” Network Description: 

The COMPUTER M system can be connected to t n r e o different systems of the 
UifiCL network. This connection allow a COMPuTERM user to work interactively 
with c t L o nil transfer files to C 0 M ? U T C R M . 


Ttui C0y<PUT£Rr-1 has one of its UZRT's connected to a 300 bps modem that 
links this micro computer tc the U t C ~ £ C L Micom port selector. This selector 
routes connection to one of the nine systems available at this network. 
Currently/ the lab has access tc the following systems: 

a CCL a : POP/10 computer running under T0PS/20 operating system. 

It has two accounts available for lab use: 
account: sosci password: ask instructor 

account: kurdahi password; ask instructor 

* GUM6Y : PQP/10 computer running under TOPS/20 operating system, 

account: smyth password: EES60L 

* VIRGIL : IBM 2/0 series 4341 running under CMS operating system, 

account: smyth password: SE56D 


rh. 

U S C - c C L : 


a i k 

* COMPUTER ;•] x 

* x 
« X x XX* X * k x * X X 


following diagram describes the interconnection between the lab and 


********** 

+ * ECLA * 

* POP/ 10 * 
********** 


********* 
.* x 

C « > * mod cm * 

k k 

k X irk XX k x X 


telephone 
<-// iane //=> 


*★ k * ★ ★ k X ★ ★ X tc 

a Micom *-+ 

* Port *~- 

* selector *- + 
************ 


************ 
■* VIRGIL * 
* IBM/ 43^1 * 
************ 


********** 

+ * GUMBY a 

A PDP/10 * 
********** 


D-i 


I 
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- C PM communication utilities. 


* TALK : Program talk.com. A1 io in t e rac t i v e communication with any 
of the £CL“sys terns mentioned above. 

* TCPCPM : Program topcpm.com. Allows text fila transfer from 5CLA or 
CUM3Y (TOPS/iG) to COMPUTERM. 

* CMSCPM : Program c.uscpm.com. Allows file transfer from 

VIRGIL (CMS) to COMPUTER M . Pile transfer is limited to S-rucord 
files only. S-record is a file format created by Motorola and 
used to allow easy inter - co:nput?r transportation of object code 
for their microprocessors. 




ill.- How to establish connection: 

1.- Boy t COMPUTERM system and turn on modem on auto - a newer mode. 

1.- Dial to SCL: 743-5Q3G or 743-7646. If you are cn campus omit dialing 
the prefix 743. 

o.- Invoke the program TALK.COM/' as soon as the carrier is detected by 
the modem: 


A > T A L K 

Wait for five 
be prompted: 


seconds and then hit <CR>. 

pvatS 



The following message will 


JSC-ECL Micoiti pert selector. 

Which system? 

3.- At thi . point/ you should select your system typing I6M/ GUM3Y 
or EiLA. 

6.- System identification will be prompted indicating that you can now 

logon to t no account. The ECL-TuPS/20 and ECL-I9M/CMS manuals provide 
a description of the main commands available in these operating 
systems. 


7.- At any point you can go back to CPM by typing control-tilde (ctrl-~). 
This w-.ll not disconnect you from ECL. Communication can be 
reestablished by invoking 7ALK.COM. 


D*-2 


3 


. - w' h i n you or o ci on > / logoff/ return to C PM ( Ctrl ) an cl turn modem off 
IV.” rile transfer to COMPUTE KM: 

So far to the documents: 

* C ommunic stion a&tuissn COMPUTERM and ISM 370/4341 ( ECL-V IR GIL) • 

* Communication oetuioan COMPUTERM and POP/10 (ECLA or GUMBY). 
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Communication o o t w e c n 


CCMPUTrRM and! I3M 370/4741 CrfCL - VIRGIL) : 


* This process requires tuio orocrairs: TALK.COM and CMSCPM.COM 

* Tc sat the connection between the two systems follow the 
following steps: 

1. ~ Plug and turn modem on and dial 5030 or 7646. 

2. - Cn CPM: execute TALK.COM 

A> T ALK 

this will cause a virtual attachment of your terminal 
to the MICCM port selector. When you select ISM and 
you are asked about your terminal type/ you should 
answer DM1320. If you want to return 'to CP My type 
Ctrl-”. This j-tui-Hr Qot log; you off nor detach you from 
the system. Communica ti'cn'-can be set again by 
repeating step 2 . 

3. - When you are don j/ log off from VIRGIL and return to 

CPM using ctrl - ”. Unplug the modem. 

: * rile transfer is limited only to one direction/ VIRGIL - > COMPUTERM 
and can only be used with 3-format record files with less than 
•» S K bytes. The /tops that must be followed -are: 

1. - jet communication between the two system using the 

•j r o c e ci u r e described above. 

2. ~ When you are ready to transfer the file/- return to CPM 

and execute TCPCPM.COM: 

A>T0PCPM XX ds s tinat icn_ f i le _in_C PM 
4>T0?CPM XX 0: TEST .03 J 

This will reconnect the two systems in the same way as 
TALK.COM does. If by some reason you have to return to 
CPM before step 3 / do it in the usual manner (ctrl-"). 

3 . - On CMS: type the CLEAR key and execute t'ne EXEC file 

CMSCPM EXEC: 
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1 

CM5CPM source filename in CMS 

" ! 

C M 5 C ? M TEST i 

ths file ext an j ion C3JCT is assumed/ ie/ TEST stands for 
TEST 0 ft J C T . The file mill b> typed to the screen as it 
gets copied. The process is slow because every time 
that the screen Is filled/ VIRGIL ujaits 40 seconds before 
continuing the transfer. At the end a message will 
b e displayed and control mill be passed to C P M . If 
you notice that the program shows abnormal behavior 
or you want to abort it/ you should 

press the RESET button in the COMPUTERS, box. Then reboot 
CPK and use TALK.COM tc retry or logoff. 

Use TALK.COM to continue or tc logoff from VIRGIL. 

If your CMS file exceeds 43K bytes. You should split 
your file in a number of files with a permissible size 
for transfer. To reconstruct the file on CPM use the 
PIP. COM utility: 

A>P I P 0ig_File = Small_File_1/.. . ./SmsllFile_n 

A>P I P A: TcST.0 2J = 3:TESTl . 0 3 J/ TE ST 2 . 03 J/ A : T EST3 . 03 J 
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APPENDIX E 


LISTING OF THE OBS SOURCE PROGRAM 


00! 3« 

o o ^ o 

03 5 3 3 
G09 03 
005 03 
0 0 5 3 3 
0073 3 
G3R03 
00933 
01333 
01133 
01233 
01303 
01903 
01533 
01509 
01703 

cmi 

019 0 3 
02030 
02133 
02233 
0233 3 
0 2933 
02530 
0 2539 
0270 3 
02333 
0293 3 
0 3303 
0513) 
03203 
0333 3 
03933 
03503 
0 3533 
03703 


P5)GU1 $IM3DJ r ( IM’ur OUTPUT); 


cosisr 

r 3 m > r i =i. ? ; 

= 13.75 ; 
G2 5 \l - 10.3? 

/: .):r7r = 2)3.3; 
Tt /E-DEITY = 2 30 .0 


:d vjstp 

: )'45T5 

run = 

T V U2 = 
r VJ5 = 
T- I M ft _ 
P H I . I M 
\l = 5? 

:* n = 

G J I 3 3 N 


= 1 . 0 ? 
= ? . o ; 
0 . 93 ; 

0-255 

0 . 53 ; 

= 23.) 
= 13.3 

57.3 ; 

= 1,0 


iplahbd; = 9 


.3 


VAR 


FU\j;-:»ETA»ltl£-l!D»U'iPOft t 5;i'U3:»3A‘H4H*GAMMAr : PEA_; 

■" n , xh , fH , xot , Y)r . y)H , ydh ♦ deltaxt , deltayt 


<r . . . 

D I L f A K H 
SI 3“IA3 


« 

, 3E. 
: 9 E A 


r h 


NEWS IGH; , 3.35I3MA 


7A2H5L-: vanes sr^inG jih r re"r to the ta 

r I =* 3 f Tiz\ 3< ♦ T.EiDXD » T) EL r , T3ELT3 , TPRI , 
f * SI , , r = SI3D ,TF5i: : 2 : VL ; 

u: v auasl: n'ns reelative to the: H3s r rune *> 


REAL 


TO THE TARGET *> 
TPRI3 


REAL 


:*l 


.. * ^ ' » J u. ^ UH I A* «. I w « 1 . » ^ ^ T U. - W 

E 3 3 » LE£3X?.EA 3X3»3~_T tDE-TD» 3 4l » ID» PSIfPSIDt PSIDOtTl REAL? 

) : ARRAY 1 1 » .13 1 0 *’ real; 

ARRAY II. .13,1.. 23 0 r REAL! 

ARRAY C 1 . • 1 0 .3 3 F REAL ? 


▲ 

or 

T* 


0393 ) 

Miyri ; hteser; 

09003 

< » j » k : integer; 
d: . r a r : r e a . ; 

09100 

0923 ) 

?si: ; 2 e a . ; 

09303 
0 9 9 0 3 
09533 

Psir : rea- ; 

09303 

09703 

B E J 1 4 (fr MAIN =*R DSR A 1 

09333 

t ;= 3 . 3 ; 

09900 

. •: ia 3 x 3 : = 3 . 3 ; 

053 3 3 

t.eoxd := 3.3 ; 

05103 

) E l r o := 3 . 3 ; 

05203 

T3E.T3 != 3.3 ; 

055 3 3 

3 H i ) := 0 . 3 ; 

C59G3 

r 3 n 1 d ;= g . o ; 

0 550 0 

3 si33 ;= 0 , 0 ; 

0 553 3 

r 5 si 3 3 ;= ).) ; 

05700 

3 s 1 3 :=o»); 

0 53 D 0 

r 3 3 1 3 ;= 0.0 ; 

0590 3 

printi ;= 3 ; 

06000 

LEA3X l- 0,0; 

OS139 

T. E A 3 < ;= 0.3 ; 

03203 

□ E L r := 3.35 
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0 6301 

csn ) 
05500 
C 5503 
C 6 7 3 1 
057 00 
0 6^31 
0733 ] 
07133 
07135 
07110 
C 7 1 1 5 
0753 3 
0 74 0 0 
0 S 3 3 3 
0 B 1 3 3 
09703 
0 9 3 3 3 
Cf 400 
095 00 
0 9 7 0 3 
039 00 

0 3 9 0 3 
09333 
0913 3 

rr\ 0 9233 
09533 

M 09409 
09533 
G 7 S 33 
09700 
093 3 3 
09903 
10 0 QO 
13133 
10? 03 
10500 

1 0 '4 0 3 

10500 
13 5 3 3 
10700 
103 00 
13930 
11000 
11133 
11233 
11303 
11»33 
11500 
11500 
11703 
119 03 
11903 
12003 
l?l 03 
122 3 3 
12300 
12433 
12405 
124 1 ja 
1250 3 ® 
12533 w 


n-.r := n, o ; 

hi := 3 . 3 : 

TI hi : = 0 » ? ; 

D : wT AT := 3.35? 

=> 3 1 r := ?73.) 7 : 1 * 3 ; 

5 s : : = 90 . 379 9 13 ; 

B I 3 * 9 1 3 := 3.3 ; (' .Hi 0- >13 U k ) 

xr := 2003 . 3 ; 
y r ;= 5 C 00 . 0 ; 

ii-jsisma ; = 4?:rxN( yt/xd-psi? 

3i ^ >4 ft h := ?o. 

3MHI := 273 . 0 7 :»h); 

wiiteln (* 04 joard aii n m mEsiEPriON shjlatih* ) ; 

W 7 HELf ; 

<3 := 0.0 ; 

yh := 0.0 ; (» imitial position 3 - he host is o»o *> 
cu :dipjm*i) m ) r he initia. .im: 3 - sight *> 
j^iTEL'J e m r :u L uz 3 - jiht ========> •♦mewsisma) ; 

1.: • mjl at 1 cm 5.0 sjidajce loop j 

W1I.E (T <= FFIklA.) DC 

3 : SIM 

<► IMITIA. t z 1 T I 3\J WI 7 HI M 7 H I ' L 3 0 p *> 

x o t ;= xr ; 

ror := vr ; 

X3H ;= <h : 

yoh := Y3 ; 

oiosn^; := jewsigma ; 


t* G'JI 31 M:E ,n 4 calculation *) 


7 4 N GE := 537 T ( S 07 ( *DT • < 3 1 ) » 3 Q 7! ’YD T- Y3H> ) ; 

ETA := SAM^A 1 -*. NEJSIGMA; 
a := g a *h a r • nensigma; 

7E-PD := (1.3 • ( ( TVELS:irY-233Cft) >7 ( VELOCI TY*COSCETA ) 3 ) ) ? 

_ a i boa := :a 3 lamjda 4 ?e.»d ; 
smmad: ;= lamso; * sismaj 
3 3 1 C := 3A7MA33 4 GUI3GN J 

(» SET F Hi OEPiVATIVE E3JATI3M3 *> 

i 3 3 := =>31 :* lO'i sr 4 ♦ d$i-:) i sr5r?si d ; 

rE = s := r 3 n: • :o.mst4*t 3 oi • :om3T5*tpsio ; h target »> 
r.EAoxo := ( r e 5 s* c om s ti • < t .:ax)7tau2 ; <* ta7get *> 

- E A CKO := (£>S*C)-:S71 - .£ A)X) 7TAJ2: 

roiLT := r a d< *■ tleaoxo * hji ; 

3 E - T := -I A 3 < 4 .E ADYD^TA J 1 S 
3 H I D := DELHC3MjT2; 

r»HiD := r 3 E . t * cgmst2 ; <* tai&-:t »> 

(* S-U M r .HUT T7E AI*C7A-T >D.w ANG.E ►> 

l- (3^140-110 > 0.0) AMO (ABMHI) > PHILIM) THEN MID := 0.0? 


H. 5 HJMT LI 1 1 r r HE TARGET AIUMfT 
IFtTMIJ^Ml > 3 . 0 ) AMD ( A 3 3 ( T => 
3 S I DO 1= { (G7A/7/ELDCI7Y) * Mt » 3 3 1 
r p s 1 33 := (<sp;.v 7 r v:L 3 cirv)'T 3 -K- 


*> 



> PHI.IM) THEM TPHI 3 
TA U 3 J 

I 0 J 7 TAJ 5 5 


= 0.3 ; 


4 ) 


l 

• 
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12 7 3 


;[r: = ) p ts 3r:o jid^^ry o 

• 

1 2^3 3 

( » 


1 293 3 
1 53 3 3 

0: -T&XT 

= Ti/:.n:riY * * c 3 sc => s i ' r > ; 


1 33 G 5 

3 E _ T A Y T 

= nz.o: ity *• d r l r ,;r » sr\i<ps T n; 


13233 

< T 

r k)t * d e l r a v r ; 


1330) 

rr 

= Y 3 T «■ deltay- ; 


13)3) 

3E.TAXM 

= »/ : l :< rt ty ► jz.i ► :)3 < 3 si >; 


L 3) 3 5 

3ELTAYH 

= VE13CITY o 5:.rA r * SI4<=SI>; 


13S0) 

< 3 

= X)H - DELTAXM ; 


1 373 3 

f H 

: Y ) H «■ CCLTAYM 5 


13233 

G A ‘SHAH 

= METANt (YH-YIH) /HH-XD-0) ; 


13303 

5 A >1 MAT 

: ARCr U< (Y7-Y3H 7t <T *X3 r >) *, 


1)333 

\1E JZi 3M4 

:= V p: TUN (( YT' YH) /( xr- YH) ) ; 


1)13 3 

SI5MA3 

:: t NI MS I GMA-OLDSIi HU 'DELTAT * 


1)20) 

1)503 

1))G3 

t* ASSI3 

4 * ) 


1)533 

i C l ft 3 : 

= _E\D< 5 

* " 

1)533 

TICl.l 3 

= r.HD< ; 


1)703 

i L l ,2 3 : 

LEVDxn ; 

1 - * .* 

1)33 3 

r i c i , 2 1 

: T.E4D<D ; 


1)903 

I C 2 »l 3 : 

3E-T ; 

0.0 

15003 

ri: 2 ,n 

= T 3 E - 7 5 

■n 2 

15133 

I C 2 ♦ 2 ) : 

):.T3 ; 

*XJ O 

^ 15203 

' 1533) 

ri:2,2i 

: T) E > T) ; 

o S 

1 : 5,13 : 

3 nr ; 

o| 

15)33 

n:3,i3 

= t=h: ; 

30 r* 

^ 15 3 0 3 

r c 3 t 2 i : 

3 hi3 ; 

O -O 

1 55 3 3 

n:3,23 

= t 3 hi o ; 

c 7? 

15703 

r o ,i 3 : 

3 s r d ; 

j> o 

15503 

no.ii 

= nsr o ? . 

r" r*i 

1593 ) 

10.23 : 

=>5133 J 

3 & 

15000 

TIC), 23 

= nsio) ; - > „ 

13133 

I [ 5 , 1 3 : 

3 3I ; 


15203 

7IC5.13 

= r =- sr ; 


15303 

IC5.23 : 

3 sr 3 ; 


15)33 

ri:5,23 

= nsro ; 


15533 

15503 

o nr:; 

0 \ r r 3M »i 


157 3 3 

j := 

l r) )0 


15333 

CCJ3: 

= i: j»i 3 «• ic j, 23- arirui 


15903 

ore J3 

: = ri: j, 13 ♦ ri: j,2 joelt at ; 


1733) 

17133 

[ * R E A 5 a 

iinment o 


1720) 

-ODK : = 

) :n; 


1 733 ) 

r l odx : 

= Drill ; 


17)00 

3 h i := o 

13 3; 


17503 

rpHi := 

3 r : 5 1 ; 


1753) 

psid := 

o: ) 3 ; 


177 00 

nsiD : = 

) re ) 3 ; 


173)3 

3 si := ) 

: 53 ; 


17933 

tp:i := 

or :53 ; 


1 ? 0 0 D 
13133 
1 3 ? 0 3 

t * p5i\rr 

DJT T 9 E PARAMETERS EVERY (3 ELTAT* 25 ) INTEGRATIONS »> 


13303 

3 rinti : 

= 3 Rt'u: ♦ i; 


13)33 

I" PRINT 

: = l THIN 


19500 

19513 

3E3IN 

writ: 

„\M»mG:r x.y.^i*. xt.yt.psit crad; 


1952) 

WRIT r 

LN{» HOST X » Y , 3 5 I * » X M , Y M » P3I-CRAD) ♦ 

CRA3, D 3IOERA0>; 

1 35 5) 

writ; 

.N(M3iT >ISMA f SHMU, 3 ;n3M1AN3»tN:ySIGMAnRA3,34MMAH. 


1 95 4 0 

WRITiL'J (*rRiC<IN3 IN- 

) 

IBS 3 3 

wr 1 r : i m ; 


13733 

F3R X t = l TO N 00 


13795 

3 I 3 1 Nl 


1 B 3 3 ) 

R I 7 1 L N < T » ICK, 

1 1 

13909 

wu r: lj { targ: 

r 

13905 

: \n; 


1 93 9 ) 

W R I T£ L M ; 


1=51 3 3 

print i: = 3 r i \ r i - 25; 


19203 



1 9 53 ) 

r := r «■ 3‘LrAr; 


19409 

:a ); 


1953 ) 

4Riri_N(«sJ:::3>»> ; 


19593 

IN). 



R.fl 3 •: t A* IT A ,MLR» 3AM0C »»R4'IS::.A,irA>CRA)»liLRDtG.i'1MA)::*CRA3> 


»::x»?),)cx3>; 

fijumitirs* *Tic<*iD*Ti;x* 23 » 3 rcK 3 > ; 


rr\ 

\ 


O o 
7 1 x 
5 c 5 

si 

It 

r r-i 

5 £5 





